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ABSTRACT 


In  Che  course  of  che  SAMSON  (Sources  of  Anbienc  Micro-Seismic ^Noise) 
Experiment  conducted  by  researchers  from  Scripps  Institution  of  Oceanography 
at  the  Coastal  Engineering  Research  Center  (CERC)  Field  Research  Facility 
(FRF)  near  Duck,  NC,  a  two-dimensional  array  of  24  bottom-mounted  pressure— ^ 
sensors  was  deployed  near  the  13 -m  depth  contour  to  detect  the  ^attenuated 
pressure  signal  arising  from  the  interaction  of  directionally  opposed  wave 
Crains . 


The  geometry  of  the  array  was  ideal  for  wind  wave  directional  spectral 
estimation  in  the  0.05  to  0.24  Hz  frequency  band.  The  formal  SAMSON  experi¬ 
ment  took  place  in  October  and  early  November  1990,  but  data  collection 
continued  until  June  1991.  Data  from  this  post -experiment  period  were  made 
available  to  CERC  and  were  subsequently  processed  with  a  360 -degree  Maximum 
Likelihood  Estimation  algorithm  to  determine  the  general  character  of  the 
local  wave  field,  including  specifically  any  indications  of  significant 
seaward-propagating  wind  wave  energy. 

This  report  contains  a  description  of  the  SAMSON  array,  the  type  of  data 
collected,  and  the  analysis  algorithm.  It  contains  displays,  both  graphic  and 
tabular,  of  relevant  wind  wave  parameters  deduced  from  the  analysis.  A  means 
of  acquiring  the  processed  data,  available  to  the  general  public  after  June 
1994,  is  also  described. 


Key  Words:  Frequency-direction  spectra 

Maximum  Likelihood  Estimation 
Wave  climate 
Wind  waves 


Jmc  Quality  inspected  2 


Accession  For  | 

ms  gra&i 

W  ! 

DTIC  TAB 

□ 

Unannounced 

□ 

Justification — 

By - 

Distribution/ 


Availability  Codej»_ 
Avail  and/or 
Special 


Diets 

ft'l 


Preface 

This  paper  provides  a  summary  of  a  aeries  of  wind-wave  frequency- 
direction  spectral  observations  made  with  a  two- dimensional,  high- resolution 
directional  wave  gage  of  opportunity,  being  the  result  of  data  shared  by 
researchers  from  Scripps  Institution  of  Oceanography  (SIO)  in  execution  of  the 
SAMSON  experiment.  The  work  vas  motivated  by  an  opportunity  to  conduct  a  360- 
degree  analysis  of  data  from  a  24 -gage  array  and  thus  find -an  encraegdinary 
obsejaaSigB  nf  dig— finally  disttttwtsd  amm  intray  Co  augment  under¬ 
standing  and  modeling  of  nearshore  processes  which  affect  coastal  engineering 
projects.  This  effort  vas  authorized  by  Headquarters,  US  Army  Corps  of 
Engineers  (HQUSACE) ,  under  Civil  Works  Research  Work  Unit  32484,  Direc¬ 
tionality  of  Waves  in  Shallow  Water,  Coastal  Flooding  Program.  Funds  were 
provided  through  the  Coastal  Engineering  Research  Center  (CERC) ,  US  Army 
Engineer  Waterways  Experiment  Station  (WES),  under  the  program  management  of 
Mr.  Charles  C.  Calhoun,  CERC.  Messrs.  John  H.  Lockhart,  Jr.;  John  G.  Housley; 
Robert  H.  Campbell;  and  James  E.  Crews  were  HQUSACE  Technical  Monitors. 

This  data  summary  was  prepared  by  Dr.  Charles  E.  Long  at  CERC's  Field 
Research  Facility  (FRF)  in  Duck,  NC,  under  the  direct  supervision  of 
Mr.  William  A.  Birkemeier,  Chief,  FRF,  and  Mr.  Thomas  W.  Richardson,  Chief, 
Engineering  Development  Division,  CERC;  and  under  the  general  supervision  of 
Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Chief  and  Assistant 
Chief,  CERC,  respectively. 

The  SAMSON  experiment,  including  the  directional  array,  was  planned  and 
executed  by  Drs.  Robert  T.  Guza  and  Thomas  H.  C.  Herbers  of  SIO.  Data 
processingjsofcware  was  written  by  Dr.  Joan  M.  Oltman-Shay  of  Quetft  Integrat¬ 
ed,  Inc.,  Seattle,  WA.  The  concept  of  continued  data  collection  after  the 
primary  SAMSON  experiment  is  credited  to  Dr.  C.  Linwood  Vincent,  CERC.  This 
document  was  edited  by  Ms.  Lee  Byrne,  Information  Technology  Laboratory,  WES. 

Commander  and  Director  of  WES  during  the  publication  of  this  report  was 
COL  Larry  B.  Fulton,  EN.  Dr.  Robert  W.  Whalin  was  Technical  Director. 
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1.  During  October  and  early  Novenber  1990,  an  experiment  involving  the 
|  detection  of  interacting  wave  trains  was  conducted  at  the  Coastal  Engineering 

^[n jfi  Research  Center  (CERC)  Field  Research  Facility  (FRF).  Known  as  SAMSON 

(Sources  of  Ambient  Micro-Seismic  Ocean  Noise),  the  experiment  was  designed  to 
detect  relatively  unattenuated  near-bottom  pressure  signals  arising  from 
nonlinear  interactions  of  wind  wave  trains  traveling  in  opposing  directions. 
These  small  but  detectable  pressure  signals  are  believed  to  initiate  ocean 
bottom  seismic  waves  that  can  propagate  long  distances  and  appear  frequently 
in  seismograms  worldwide.  In  detecting  both  presence  of  interacting  wave 
trains  and  resulting  pressure  signals ,  results  of  the  SAMSON  experiment  can  be 
used  to  test  hypotheses  concerning  the  initiation  of  such  seismic  waves  near  a 
coast  where  nearshore  bathymetry  can  act  as  a  partial  reflector  of  ocean  wind 
waves  and  thereby  induce  the  requisite  opposing  wave  trains.  The  experiment 
was  designed  and  executed  under  the  direction  of  Drs.  Robert  T.  Guza  and 
Thomas  H.  C.  Herbers  of  the  Center  for  Coastal  Studies  (CCS)  at  Scripps 
Institution  of  Oceanography  (SIO) ,  University  of  California,  San  Diego  (UCSD) , 
with  experiment  site,  staging  area  and  logistical  support  provided  by  the  FRF. 

2.  The  primary  instrument  in  the  SAMSON  experiment  was  a  two-dimension¬ 
al  array  of  24  bottom-mounted  pressure  gages  located  in  a  nominal  water  depth 
of  13  m.  By  design,  this  depth  is  deeper  than  is  normally  used  for  measure¬ 
ment  of  ocean  wind  waves.  The  idea  was  to  let  the  water  column  act  as  a 
filter  for  pressure  signals  arising  directly  from  wind  waves  above  a  certain 
cut-off  frequency.  In  that  way,  any  detectable  pressure  signal  above  the  cut¬ 
off  frequency  could  be  interpreted  as  being  from  the  wave-wave  interactions 
sought  in  the  experiment.  At  frequencies  below  the  cut-off,  the  pressure 
gages  could  readily  detect  wind  wave  signals.  As  a  result,  the  array  of  gages 
was  an  excellent  directional  detector  for  waves  in  the  lower  part  of  the  wind 
wave  frequency  band,  from  about  0.05  to  0.24  Hz. 

3.  The  SAMSON  pressure  gage  array  was  deployed  in  September  1990  and 
the  formal  experiment  took  place  in  October  and  early  November.  Scheduling 
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recovery  of  Che  gages  immediately  after  the  experiment  vas  deemed  usvise 
because  of  Che  onset  of  winter  and  the  high  likelihood  of  poor  working 
conditions.  Hence,  the  gages  were  left  in  place  uncil  Che  following  spring 
and,  because  it  was  a  fully  functioning  instrument  sec,  data  were  collected 
for  as  long  as  the  gages  continued  Co  funcclon,  chough  this  was  not  expected 
to  be  for  very  long  given  Che  harshness  of  the  environment^  in  which  thay~WTe 
deploy edT  Curiously,  22  of  Che  24  gages  continued  to  function  perfectly  until 
the  first  part  of  June  1991,  at  which  time  the  gages  were  recovered.  Data 
from  this  post-experiment  phase  of  instrument  deployment,  mid-November  1990  to 
the  end  of  May  1991,  have  been  made  available  to  the  FRF  for  directional 
analysis.  This  report  addresses  the  natury^and  results  of -thee  analysis . 

4.  Because  of  the  volume  of  data  involved,  considerable  time  has  been 
saved  by  performing  a  preliminary,  fairly  rudimentary  analysis  to  obtain  the 
general  character  of  the  observed  directional  energy  distributions ,  as 
described  by  a  few  simple  parameters.  More  detailed  analysis  can  be  performed 
at  a  later  date  on  individual  cases  or  sequences  of  cases  that  are  of  most 
interest.  The  high-resolution  directional  estimation  algorithm  known  as 
Maximum  Likelihood  Estimation  (MLE)  (Davis  and  Regier,  1977)  is  used  here  for 
the  preliminary  analysis  because  it  is  rather  quick  to  compute  and  is  adequate 
to  describe  the  general  directional  nature  of  the  observations .  A  more 
refined  but  more  time-consuming  algorithm,  such  as  Iterative  Maximum  Likeli¬ 
hood  Estimation  (IMLE)  (Pavka,  1983),  can  be  used  later  to  increase  resolution 
on  select  cases. 

5.  Because  the  SAMSON  array  is  two-dimensional,  it  can  detect  wave 
energy  propagating  in  any  direction.  This  allows  estimation  of  energy 
propagating  seaward  and  avoids  the  180-degree  ambiguity  present  in  analysis  of 
data  from  the  FRF  linear  array  directional  wave  gage  reported  previously  (Long 
and  Oltman-Shay  1991,  Long  1991a,  1991b).  In  this  report,  seaward  propagating 
energy  is  called  "reflected"  energy  because  reflections  of  incident  waves  from 
nearshore  bathymetry  is  its  most  likely  cause,  though  ether  pruiasser  may 

Determination  of  reflected  directional  spectra  is  of 


inori 

possibly -ba  present . 


enormous  interest  in  coastal  research  because  of  myriad  processes  that  are 
possible  theoretically  if  significant  reflection  is  present.  Hence,  this 
report  is  the  first  to  describe  a  well- resolved  estimation  of  the  full  360- 
degree  wave  climate  at  the  FRF. 
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6.  The  beginning  text  of  (JlUr^ocumenC  describes  and  clarifies  the 


substantial  information  ritntrHnnrl  fri  appgrrtttCTs^J  A  brief  overview  is 
given  of  the  measurement  site,  instrumentation,  data  collection,  and  method  of 
directional  spectral  estimation.  These  subjects  are  described  in  other 
publications  to  which  the  reader  is  referred  for  greater  detail.  Following 
the  overview  is  a  description  of  the  archived  frequency* direction  spectra  and 
some  characterizing  bulk  parameters  that  can  be  derived  from  them.  Appendix  A 
is  a  listing  of  these  characterizing  parameters  and  is  Intended  to  be  used  as 
a  EStfPf^catalog  of  the  set  of  spectra.  Appendix  B  contains  graphs  of  time 
series  of  some  of  these  parameters  as  a  pictorial  augmentation  of  the  informa¬ 
tion  in  Appendix  A.  Appendix  C  is  a  derivation  of  the  MLE  used  in  this 
analysis.  It  is  intended  to  augment  the  somewhat  brief  derivations  normally 
found  in  the  literature,,  and -thereby  peiliaps  augment  the  reader  *s  ■tuiUei’SCahiling 
of— this  method  »f  directional'  estimation. 
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Measurement  Site 


7.  As  shown  in  Figure  1,  the  FRF  is  located  on  the  barrier  island  chain 
of  coastal  North  Carolina.  A  detailed  description  of  the  layout,  function, 

^  and  capabilities  of  the  FRF  is  given  by  Birkemeier  et  al.  (1985).  A  detailed, 
quantitative  description  of  the  climate  at  the  FRF,  as  determined  from  its 
arsenal  of  instrumentation  and  for  the  ten  years  of  Its  existence,  is  given  by 
Leffler  et  al.  (1990).  Of  particular  relevance  to  directional  wave  studies 
are  the  wave -steering  bathymetry  and  wave -generating  winds. 

8.  As  regards  the  former,  the  coastline  in  the  vicinity  of  the  FRF  is 

nearly  straight  for  several  tens  of  kilometers  north  and  south  (Figure  1) .  It 
is  oriented  such  that  a  shore-normal  line  (directed  seaward)  is  very  nearly 
70  deg  from  true  north.  Waves  and  onshore  winds  can  approach  this  site  along 
an  easterly  180-deg  arc  from  340  to  160  deg  true.  The  adjacent  continental 
shelf  is  wide,  relatively  shallow,  and  of  somewhat  complex  bathymetry.  The 
direction  of  nearest  approach  of  the  100-m  Isobath,  which  indicates  the  shelf 
break,  is  10  to  15  deg  south  of  east  and  is  about  80  km  distant.  A  typical 
bottom  slope  for  the  shelf  is  1  m/km,  but  this  is  interrupted  by  numerous 
features  of  1-  to  10-km  horizontal  scales  and  10-m  vertical  scales  scattered 
irregularly  across  the  shelf.  >  -r~ 
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9.  Within  a  few  kilometers  of  the  FRF,  the  offshore  bathymetry  is  more 
regulagf  with  isobaths  nearly  shore -parallel  and  a  bottom  slope  of  about 

2  m/km  (Figure  2).  Some  irregularities  exist.  Within  about  300  m  of  the 
shore,  there  exists  a  complex  and  mobile  bar  system  (Birkemeier  1984).  Waves 
and  currents  have  created  some  irregular  bathymetry  in  the  vicinity  of  the  FRF 
research  pier  which  extends  about  500  m  offshore  (Miller,  Birkemeier,  and 
DeWall  1983). 

10.  The  site  is  subject  to  a  variety  of  climates  which  gives  rise  to  a 
diverse  set  of  directional  wave  conditions.  Primary  sources  of  high-energy 
waves  4re  winds  associated  with  hurricanes  and  frontal  passages.  Though 
Hurricane  Lily  moved  up  the  middle  Atlantic  Ocean  during  the  formal  SAMSON 
experiijent,  no  other  hurricanes  affected  the  FRF  wave  climate  in  the  period 
covered  by  this  report.  Low-pressure  weather  fronts,  of  which  several  crossed 


Figure  2.  Approximate  nearshore  bathymetry  and  SAMSON  array 
location  (seawardmost  sec  of  dots) 
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Che  FRF  site  during  this  reporting  period,  were  typically  oriented  northeast- 
southwest,  with  strong  wave -generating  winds  coming  from  the  northeast. 

11.  As  a  consequence  of  the  complex  offshore  bathymetry  and  the  close 
proximity  of  wave -generating  wind  systems,  it  is  expected  that  directional 
wave  spectra  will  exhibit  a  rather  broad,  complex  structure.  This  expectation 
has  been  realized  in  studies  at  the  FRF  using  a  linear  array  directional  wave 
gage  in  8  m  water  depth  (Long  and  Oltman-Shay,  1991).  Excepting  well-sorted 
waves  from  distant  storms  and  waves  in  the  Initial  growth  stages  under  strong 
steady  local  winds,  most  directional  spectra  were  found  to  be  rather  broad. 
Because  of  the  180-degree  ambiguity  inherent  in  a  linear  array,  it  is  possible 
that  part  of  this  broadness  could  be  due  to  reflected  wave  energy,  which  is 
folded  over  to  appear  as  incident  wave  energy  in  the  linear  array  data 
processing  algorithm.  The  ability  to  perform  a  360-degree  wave  analysis  with 
the  SAMSON  array  will  help  to  ascertain  the  extent  of  the  error  induced  by  the 
180-degree  ambiguity  of  the  linear  array,  in  addition  to  estimating  the  amount 
of  seaward  propagating  energy. 


Instrumentation 

12.  The  primary  instrument  in  this  study  was  a  high-resolution,  two- 
dimensional  array  directional  wave  gage.  It  consisted  of  two  parts.  The 
first  was  an  array  of  sensors  that  sampled  sea- surface  displacement  at  several 
points  in  (horizontal)  space.  The  second,  described  in  the  following  section 
on  data  processing  and  in  Appendix  C,  Is  the  mathematical  treatment  of  these 
data  to  obtain  estimates  of  wave  directionality. 

Basic  array 

13.  The  SAMSON  array  consisted  of  24  pressure  gages  buried  about  0.1  m 
below  the  bottoa  in  the  vicinity  of  the  13 -a  isobath  about  1700  m  offshore  and 
slightly  to  the  north  of  the  research  pier.  This  is  shown  in  Figure  2  as  the 
seawardmost  set  of  dots.  In-additinn  rn  the  design  requirements  mewtiorrgcf  ~ 
above ffi;s  location  satisfies  three  constraints.  First,  it  is  generally 
outside  the  surf  zone  so  that  linear  wave  theory  is  applicable  in  data 
processing.  Second,  it  is  in  water  shallow  enough  that  signals  from  4-sec 
waves  are  detectable  above  background  noise  at  the  bottom-mounted  gages. 

Third,  it  is  away  from  the  Irregular  isobaths  around  the  pier  and  in  the 


8 


nearshore  bar  ay seen,  which  helps  ainiaize  bathyaecrleally  Induced  inhomo gene 
lties  in  Che  wave  fluid. 

14.  Figure  j  is  an  enlarged  view  of  Che  array  layouc  and  shows  gage 
spacing  as  well  as  the  gage  numbering  scheme.  Spacing  beeveen  Che  gages  in 
Che  array  corresponds  reasonably  well  to  the  array-design  criterion  posed  by 
Davis  and  Regier  (1977)  chac  every  gage  pair  have  a  unique  separation  or 
.patial  lag.  There  are  a  large  nuaber  of  redundant  lags  built  into  the  array 
so  chat  multiple  estimates  of  cross -spectral  densities  can  be  made  to  check 


SAMSON  13-m  Array 


-128  -»« 


-32  0  32 

Longshore  (m) 


98  128 


Figure  3.  Layout  and  numbering  scheae  for  SAMSON  array 
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for  spatial  inhomogeneity  in  the  wave  field  and  to  reduce  system  noise.  There 
are  a  much  greater  number  of  unique  lags  with  which  to  fill  in  the  spatial 
cross -spectral  matrix  from  which  wave  direction  estimates  are  made.  Minimum 
gage  spacing  is  4  m,  maximum  spacing  (the  greatest  dimension  of  the  array)  is 
256  m,  and  intermediate  gage  spacings  are  in  multiples  of  4  m.  Unfortunately, 
gages  12  and  14  were  inoperative,  and  the  loss  of  gage  12  reduced  the  long¬ 
shore  dimension  of  the  array  to  192  m,  at  some  cost  to  direction  resolving 
ability  for  the  longest  (0.05  Hz)  waves  measured.  It  is  noted  that  gage  19 
appears  to  be  offline  from  the  basic  gage  pattern,  and  it  was,  in  fact,  jkJ'- 

slightly  misplaced.  However,  this  does  not  make  it  ineffective.  In  ML E  V- 

processing,  it  is  only  necessary  for  the  true  gage  positions  to  be  known,  and 
not  that  they  conform  to  any  particular  pattern.  The  SAMSON  gage  positions 
were  surveyed  to  an  accuracy  of  better  than  1  m  after  the  gages  were  in  place. 

At  worst,  this  uncertainty  would  result  in  an  error  of  about  1.8  degrees  for 
the  shortest  (0.24  Hz)  waves  considered  and  about  0.6  deg  for  the  longest 
(0.05  Hz). 

Sub -arrays 

15.  At  sufficiently  large  lags,  signals  from  a  natural  wave  field  tend 
to  become  uncorrelated.  As  a  result,  cross-spectra  between  widely  spaced 
gages  tend  to  become  mostly  noise.  At  very  small  lags,  correlations  and 
cross -spectra  are  difficult  to  discern  from  those  at  zero  lag  and  so  tend  to 
become  somewhat  redundant.  The  terms  large  and  small  are  ganarnilly  related  to 
the  wavelengths  of  waves  under  consideration.  To  avoid  both  redundancy  at 
small  lags  and  the  problem  of  noise  at  large  lags,  the  SAMSON  array  data  were 
processed  by  using  subsets  of  gages,  or  sub-arrays,  based  on  ranges  of 
wavelengths  being  considered. 

16.  Figure  4  illustrates  the  sub-arrays  and  the  ranges  of  wave  frequen¬ 
cy  used  with  each  sub-array  used  in  the  present  analysis.  Note  that  gages  12 
and  14,  from  which  reliable  data  were  not  obtained,  are  not  shown.  For  the 
longest  waves,  all  the  working  gages  were  used  except  those  which  introduced 
the  shortest  (4-m)  lags,  i.e.,  gages  5  and  22.  This  left  20  gages  for 
processing  data  in  the  frequency  f  range  0.05  <  f  (Hz)  <0.10  or  wave¬ 
length  X  range  103  <  X  (m)  <  220  (assuming  a  nominal  water  depth  of  13  m) , 
and  Included  the  longest  possible  lags. 

17.  For  shorter  wavelengths,  the  selection  criterion  was  that  the 
longest  dimension  of  the  sub-array  not  exceed  about  three  wavelengths  of  the 
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Figure  4.  Sub-arrays  and  their  associated  frequency  bands 


shortest  waves  in  a  given  band.  For  the  low  aid- range  of  frequencies, 

0.X0  <  f  (Hz)  <0.16  ,  corresponding  to  wavelengths  55  <  A  (a)  <  103  ,  this 
procedure  eliminated  gage  19  from  analysis,  leaving  19  gages  in  the  sub-array. 
For  the  high  aid-range  of  frequencies  0.16  <  f  (Hz)  <0.20  ,  or 

38  <  A  (a)  <  55  ,  gages  1,  2,  13,  and  21  vere  eliminated,  but  the  short-lag 

gages  5  and  24  vere  Included  for  a  total  of  17  gages.  For  the  highest 
frequencies,  0.20  <  f  (Hz)  <  0.24  ,  27  <  A  (m)  <  38  ,  gages  3,  11,  and  20 

were  excluded  froa  the  previous  set,  leaving  14  gages. 
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18.  Each  pressure  gage  is  a  Setra  Systems  Model  280E  capacitive- type 
transducer  having  a  range  of  0  to  50  psla.  The  manufacturer's  specifications 
indicate  an  accuracy  of  O.llt  of  full  scale  or  the  pressure  equivalent  of 
±0.037  a  of  water.  Analog  pressure  signals  were  collected  in  an  underwater 
terminal  where  they  were  discretized  to  12 -bit  binary  form,  multiplexed,  and 
encoded  for  transmission  down  a  single  data  line.  Digitization  of  the  signal 
to  12-bit  binary  form  resulted  in  a  discretization  uncertainty  of  the  equiva¬ 
lent  of  0.002  m  of  water.  At  the  landward  end  of  the  data  line,  the  encoded 
signal  was  decoded  and  sent  to  a  parallel  port  on  an  AST  286  personal  computer 
for  data  acquisition.  During  the  post-experiment  period,  data  files  from  the 
personal  computer  were  periodically  transferred  to  the  higher -capacity  hard 
disks  of  the  FRF  VAX  11/750  computer  for  accumulation  and  eventual  archiving 
on  magnetic  media. 


Data  Collection 

19.  Signals  from  each  of  the  24  pressure  gages  were  sampled  at  4  Hz, 
and  stored  in  records  of  768  samples.  A  total  of  1,280  records  or  10,240  data 
points  were  collected  for  each  gage  for  an  elapsed  time  of  2  hr  50  min  40  sec 
for  each  collection.  Starting  times  for  collections  during  the  main  SAMSON 
experiment  were  daily  at  0100,  0400,  0700,  1000,  1300,  1600,  1900,  and  2200 
Eastern  Standard  Time  (EST)  and  this  pattern  was  continued  into  the  post- 
experiment  period  until  1  December  1990.  At  that  time,  collection  start  times 
were  set  to  coincide  with  routine  FRF  observations  (Blrkemeier  et  al.  1985), 
which  occur  daily  at  0100,  0700,  1300,  and  1900  EST.  This  collection  pattern 
was  maintained  until  the  gages  vere  disconnected  for  retrieval  on  1  Jun  1991. 

20.  In  the  analysis  presented  here,  only  the  first  2  hr  16  min  32  sec 
(eight  records  of  4,096  points)  of  data  from  each  collection  was  used.  This 
procedure  resulted  in  data  sets  that  coincided  both  in  start  times  and 
duration  with  other  permanent  directional  wave  gages  at  the  FRF.  It  also 
allowed  consistency  in  processing  some  collections  that  vere  terminated  early 
to  allow  time  for  file  transfers  from  the  data  collection  computer  to  the  FRF 
VAX  11/750  hard  disks. 
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Bata  Exscgaalng 


21.  Conversion  of  oeasured  cine  series  Co  esciaaces  of  frequency- 
direction  spectra  requires  products  of  frequency  spectral  estimates  from  all 
the  various  gages  in  the  array.  For  final  results  to  be  accurate,  raw  input 
data  must  be  of  exceptionally  high  quality  so  that  spiky  or  drifty  data  from 
one  gage  do  not  contaainate  products  of  results  froa  the  other  gages.  Hence, 
Che  procedure  for  data  processing  is  to  check  the  raw  data  for  errors, 
estimate  the  frequency-direction  spectrum,  and  then  compute  some  bulk  para¬ 
meters  with  which  to  characterize  results. 

Error  checking 

22.  Because  multiple  gages  were  deployed  in  an  (assumed)  uniform  sea, 
certain  properties  of  raw  data  from  the  pressure  gages  should  be  identical.  A 
convenient  method  of  error  checking  is  to  compute  frequency  spectra  of  the  raw 
time  series  from  each  of  the  22  functioning  gages  and  overplot  the  results  on 
a  single  graph.  This  metliuUmnr  used  by  PrS  .'  Cuza  and  HBfb*ri»  during  the  main 

rirper-tmanc  nnri  haa  elan  bean  inocrp^rated  fftr  flnr  present  analysts .  Because 

of  the  depth  at  which  the  SAMSON  gages  were  deployed,  very  little  direct  wind 
wave  signal  is  expected  in  the  frequency  range  froa  about  0.3  Hz  out  to  the 
Nyquist  frequency  of  2.0  Hz.  The  frequency  spectrum  in  this  frequency  band 
should  primarily  indicate  system  noise,  which  should  be  about  the  same  for 
each  gage.  Excessively  spiky  data  froa  one  gage  will  raise  the  noise  level 
for  that  gage  relative  to  the  noise  levels  froa  the  other  gages  so  that  its 
signature  will  be  distinct  in  a  collection  of  overplotted  spectra.  Strong 
drifts  in  the  signal  can  appear  in  a  variety  of  ways  in  corresponding  spectra, 
but  will  generally  give  a  distinct  signal  at  low  frequencies.  In  the  pass 
band  of  wind  wave  frequencies  for  which  directional  estimates  are  computed 
(here  0.05  to  0.24  Hz),  one  expects  the  frequency  spectra  to  be  nearly 
Identical . 

23.  Figure  5  is  an  example  of  one  set  of  overplotted  frequency  spectra. 
Semi -logarithmic  coordinates  have  been  used  and  pressure  gage  voltages  have 
been  converted  to  equivalent  heights  of  a  static  water  column,  but  no  other 
transformations  have  been  used.  As  can  be  seen,  signals  in  the  wind  wave 
frequency  pass  band  are  all  very  nearly  alike,  indicating  that  all  gages  are 
functioning  reasonably  well.  The  noise  floor  at  high  frequencies  is  very  low 
relative  to  the  wind  wave  signal  and  is  nearly  uniform  for  all  gages.  The 


■pikes  in  the  spectre  *t  high  frequencies  ere  due  to  low  levels  of  commercial 
60  Hz  power  line  noise  present  in  the  system.  The  column  of  numbers  in  the 
right  asrgin  is  the  gage  number  associated  with  each  spectrum.  If  a  particu¬ 
lar  gage  becoaes  acre  noisy  than  the  others,  its  identifying  ttuaber  can  then 
be  readily  determined.  Because  the  signal  aean  is  removed  prior  to  spectral 
estiaation,  it  would  not  normally  appear  in  a  spectral  plot.  Signal  Beans  are 
also  useful  indicators  of  misbehaving  gages,  so  they  are  written  out  in  the 
two  rows  of  nuabers  below  the  plot  header. 

24.  Initial  tests  with  this  kind  of  checking  routine  indicated  that  the 
SAHSON  gages  and  data  collection  systea  were  extremely  reliable.  Consequent¬ 
ly,  overplotted  spectra  were  not  computed  for  every  run,  but  were  computed  for 
the  first  run  of  each  day.  At  no  tine  during  the  post -experiment  period  were 
any  of  the  22  functioning  gages  found  to  have  any  significant  errors. 
Frecuencv-direction  spectra 

25.  Estimation  of  frequency-direction  spectra  is  done  in  four  parts. 
First,  time  series  of  pressure  data  from  each  gage  are  Fourier  transformed  to 
the  frequency  domain.  Second,  these  transforms  are  converted  to  sea-surface 
displacement  transforms.  Third,  cross-spectra  of  sea -surface  displacement  are 
computed  between  all  unique  gage  pairs  for  each  frequency.  Finally,  an 
estimate  is  made  of  a  directional  distribution  of  vave  energy  that  corresponds 
to  the  computed  spatial  variation  in  cross -spectral  density  for  each  frequen¬ 
cy. 

26.  The  Fourier  transform  is  conventional.  A  32, 768 -point  tine  series 
is  divided  into  15  half -overlapping  segments  of  4,096  points.  Segments  are 
tapered  with  a  Kaiser-Bessel  window  (a  modified  Bessel  function  of  the  first 
kind,  compensated  uniformly  for  loss  of  variance  due  to  windowing)  and  fast 
Fourier  transformed.  An  intermediate -resolution  transform  is  found  by 
averaging  the  15  transformed  segments,  frequency  by  frequency.  Final  trans¬ 
forms  are  found  by  then  averaging  results  over  10  adjacent  frequency  bands. 
Final  resolution  bandwidth  is  0.00976  Hz,  and  degrees  of  freedom  are  at  least 
150  (assuming  8  contiguous  segments  and  Ignoring  any  gain  from  lapped  seg¬ 
ments).  Transform  estimates  are  retained  for  20  frequency  bands  with  band- 
center  frequency  ranging  from  0.054  to  0.240  Hz. 

27.  Conversion  of  pressure  signals  at  depth  to  water -surf ace  displace¬ 
ment  is  done  through  the  linear  wave  theory  pressure  response  factor  as 
described  in  the  Shore  Protection  Manual  (SPM)  (1984).  After  this  conversion, 
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complex  cross-spectra  in  the  fora  of  coincident  and  quadrature  spectra  are 
computed  in  the  conventional  way  (Bendat  and  Piersol  1971,  Jenkins  and  Watts 
1968)  between  all  unique  gage  pairs.  Cross -spectral  estimates  at  a  given 
frequency  are  then  stored  in  a  matrix  of  coaplex  numbers,  the  indices  of  which 
can  be  used  to  address  data  arrays  containing  gage  position  for  the  computa¬ 
tion  of  gage  separation  distance,  or  lag  space,  in  preparation  for  directional 
spectral  estimation  at  that  frequency. 

28.  In  this  preliminary  analysis  of  the  SAMSON  data,  conversion  of 
cross-spectral  patterns  in  lag  space  to  directional  spectra  is  done  with  the 
Maximum  Likelihood  Estimation  (KLE)  algorithm  described  by  Davis  and  Regier 
(1977).  The  algorithm  is  also  described  in  application  to  data  from  heave- 
pitch-roll  buoys  by  Oltman-Shay  and  Guza  (1984).  Neither  of  these  two 
references  includes  all  the  steps  in  the  derivation,  so,  for  pedantic  com¬ 
pleteness,  a  more  complete  outline  of  the  derivation  of  MLE  is  given  in 
Appendix  C. 

29.  The  algorithm  used  here  made  a  spectral  estimate  for  every  1  deg  of 
azimuthal  direction.  These  were  smoothed  over  2-deg  arcs  to  create  180  bins 
covering  a  complete  circle  with  which  to  describe  the  directional  distribution 
of  wave  energy. 

30.  /The  primary  result  of  data  processing  is  an  estimate  of  the 
discrete  frequency-dlrectionSapectrum  S(f„,0B)  ,  which  represents*  the  vari¬ 
ance  of  sea- surface  displacement  per  frequency  resolution  bandwidth  df 


(-  0.00976  Hz)  per  direction  resolution  arc  it  (-2  deg),  where  f,,  is  the 
n**  of  N  -  20  discrete  frequencies  and  §m  is  the  a*  of  M  -  180 

\ 

discrete  directions.  In  this  work,  direction  is  considered  to  be  the  angle 
from  which  wave  energy  is  coming,  measured  counterclockwif*  from  shore  normal . 
Also,  it  is  pointed  out  that  parametric  definitions  below  are  given  In  terms 
of  these  discrete  definitions  rather  than  in  terms  of  the  continuous  functions 
they  approximate. 

31.  Numerical  values  of  S(fa,#B)  can  range  over  many  orders  of 
magnitude,  depending  on  the  amount  of  energy  in  a  given  frequency  band  and 
direction  arc,  and  this  can  require  space -consuming  formats  for  archiving 
data.  To  simplify  this  problem,  frequency- direction  spectra  can  be  saved  in 


For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation  / 
(Appendix  D) .  x 

, wr 
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(1) 


Che  form  of  directional  distribution  functions  D(fn,  5.)  defined  by 


S(fa) 


where  S(fn)  is  the  frequency  spectral  density  at  frequency  .  The 
directional  distribution  function  has  units  of  deg'1,  and  its  integral  with 
respect  to  direction  over  all  directions  is  unity. 

32.  The  frequency  spectrum  in  Equation  1  represents  the  total  over  all 
directions  of  sea*  surface  variance  per  frequency  bandwidth  and  is  defined  in 
terms  of  the  frequency- direction  spectrua  by 


s(fa)  m£ 


(2) 


Note  that  this  is  the  sane  thing  as  a  conventional  frequency  spectrua  that 
would  result  from  a  tine  series  of  sea- surface  displacements  at  a  single  point 
in  space.  Because  it  is  an  iritegral^of  the  frequency -direction  spectrum,  it 
is  called  the  ^hntegrateB^ frequency  spectrum. 

33.  A  directional  analog  of  the  frequency  spectrum  is  the  Integrated 
direction  spectrum,  found  by  summing  the  frequency-direction  spectrum  over  all 
frequencies  for  a  fixed  direction  arc.  It  is  computed  from 


"  E  S(fu,$m)df 

n»l 


(3) 


Figure  6  shows  one  way  to  display  the  frequency-direction  spectrum  and  the 
corresponding  integrated  frequency  and  integrated  direction  spectra. 

Bulls  parameters 

34.  Several  parameters  have  been  computed  to  characterize  the  observed 
spectra.  Five  elementary  types  of  parameter  have  been  used  here:  water 

depth,  characteristic  wave  height,  peak  frequency  (or  its  Inverse,  peak _ 

period) ,  peak  direction,  and  reflection  coefficient^rThare  is  more  than  one 


way  to  define  some  of  these  parameters  r~  so  definitions  are  presented  here. 

35.  Water  denth.  Water  depth  is  found  from  the  average  of  the  mean 
depths  from  all  22  gages  in  the  array.  It  represents  a  water  level  that 
varies  with  low  frequency  processes  like  tides.  Because  variations  in  meanN 
water  level  adjust  intersections  of  the  vater  surface  with  nearshore  bathy-/ 
me  try,  it  may  be  an  important  variable  in  accounting  for  the  amount  of 
observed  reflected  energy.  I  ± 
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36.  Characteristic  wave  height.  Characteristic  wave  heights  from 
spectral  observations  are  most  frequently  given  as  ,  which  is  four  times 
the  standard  deviation  of  sea-surface  displacement.  Zt  can  be  determined  from 
the  volume  under  the  frequency-direction  spectrum  by  the  equation 

SL  •  16  E  E  S(fa,9m)  dBdf  (*> 

a»l  a>l 

It  can  be  found  equivalently  from  the  integrated  frequency  spectrum  by 

s 

-  16  £  S(fa)  df  (5) 

#■1 


which  is  its  conventional  definition,  or  from  the  integrated  direction 
spectrum  by 

j£«  16  E  S(em)d8  <6> 

37.  Because  the  analysis  yields  a  360-deg  distribution  of  wave  energy, 
it  is  reasonable  to  distinguish  between  Incident  wave  energy  and  outgoing  wave 
energy,  which  here  is  arbitrarily  called  reflected  wave  energy.  A  character¬ 
istic  incident  wave  height  can  be  identified  from  the  volume  under  the 

frequency-direction  spectrum  for  waves  approaching  from  directions  in  the 
range  -90  deg  <  9  <  90  deg  or,  equivalently,  the  area  under  the  integrated 
direction  spectrum  over  the  same  range  of  directions.  Discrete  directions  in 
the  present  analysis  are  defined  by 

9m  -  180  -  2(m  -  1)  (7) 

If  half  the  energy  in  the  90-deg  and  -90 -deg  bins  and  all  the  energy  in  the 
intervening  directional  bins  is  attributed  to  incident  energy,  a  characteris¬ 
tic  Incident  wave  height  is 

133 

*L.i  -  4S(S»«)  d8  *  E  s(*-> da  *  45(*ij6>  de  (8) 

2  2 

If  the  other  half  of  the  energy  in  the  90-deg  and  -90-deg  bins  and  the 
remainder  of  the  offshore  propagating  energy  is  attributed  to  reflected 
energy,  a  characteristic  reflected  wave  height  is  given  by 


19 


(9) 


45  ISO 

i£  t  -  E  s<*.> d#  ♦  rf#  ♦  E  *<'.>  di 

mml  2  2  e-U7 

38.  Peak  frequency.  Peak  frequency,  which  ha*  Che  generic  notation 

fp  ,  can  be  defined  aC  least  two  ways.  One  way  is  to  find  the  frequency  (and 
direction)  at  which  the  frequency-direction  spectrum  is  maximum.  Another  way 
is  to  find  the  frequency  at  which  the  integrated  frequency  spectrum  is 
maximum.  This  is  the  more  conventional  definition  and  is  the  basic  definition 
used  here.  As  with  characteristic  wave  height,  a  distinction  is  made  between 
peak  frequencies  associated  with  incident  and  reflected  wave  energy.  The  peak 
frequency  of  the  incident  wave  spectrum  fpit  is  defined  as  the  frequency  of 
the  maximum  of  the  integrated  incident  frequency  spectrum  S^f.)  defined 
following  the  notation  above  for  as 

133 

dt  *  E  dS  ♦  7 *(*..*«•>  d*  (10) 

2  2 

In  a  similar  manner,  the  peak  frequency  of  the  reflected  wave  spectrum  f  s 
is  defined  as  the  frequency  of  the  maximum  of  the  integrated  reflected 
frequency  spectrum  S*(fn)  given  by 

43  1*0 

st(fa)  -E  S(fa,ea)de  ♦±s(fm.i„i)  ♦  £  s(fn,8m)de  UU 

••1  2  ••137 

The  distinction  between  fpit  and  fp>c  may  be  istportant.  Unless  there  are 
very  strong  reflections,  the  incident  peak  frequency  is  virtually  the  same  as 
the  peak  frequency  one  would  obtain  from  a  time  series  at  a  single  point  in 
space.  If  the  beach  is  uniformly  reflective  (the  same  for  all  frequencies), 
reflected  peak  frequency  will  be  the  same  as  Incident  peak  frequency.  If  the 
two  frequencies  are  not  equal,  the  beach  may  be  significantly  more  reflective 
for  certain  frequencies  other  than  the  peak,  a  situation  of  no  small  dynamic 
importance  in  nearshore  processes. 

39.  Peak  period.  Peak  period  Tp  is  the  characteristic  wave  period 
associated  with  spectral  peak  frequency.  As  simply  inverses  of  peak  frequen¬ 
cies,  there  are  two  peak  periods  associated  with  a  directional  spectrum 
corresponding  to  the  two  peak  frequencies  defined  above.  The  Incident  peak 
period  is  Tp>1  -  I/fp  l  and  the  reflected  peak  period  is  Tp>r  -  l/fp>r  . 

40.  Peak  direction.  Peak  direction  is  the  direction  at  which  the  most 
energy  occurs.  Given  the  generic  symbol  tp  ,  it  is  defined  analogously  to 
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peak  frequency  through  the  Integrated  direction  spectrum.  As  with  the  other 
parameters  there  is  an  incident  peak  direction  #PiL  and  a  reflected  peak 
direction  #p>r  .  Incident  peak  direction  is  the  direction  of  the  maximum  of 
the  Integrated  incident  direction  spectrum,  which  is  Equation  3  restricted  to 
those  directions  for  which  -90  deg  <  9  <  90  deg  or  for  which  index  m  is  in 
the  range  46  <  a  <  136  .  Reflected  peak  direction  is  the  maximum  of  the 
Integrated  reflected  direction  spectrum,  again  Equation  3  but  now  restricted 
to  offshore  propagation  directions.  This  restriction  leads  to  two  ranges  (in 
the  present  definition)  of  direction  -180  deg  <  $  <  -90  deg  and  90  deg  <  8 
<  180  deg  or  of  direction  indices  1  <  a  <  46  and  136  <  m  <  180  .  These 
measures  of  peak  direction  indicate  only  the  single  direction  cells  in  each 
direction  range  that  contain  the  most  energy.  If  much  depth-  and  frequency - 
dependent  refraction  is  occurring  (as  appears  to  be  occurring  in  Figure  6, 
e.g.),  the  Integrated  direction  spectra  are  somewhat  smeared  across  the 
direction  axis.  Peak  directions  thus  found  may  not  correspond  then  with  any 
peaks  in  the  full  frequency-direction  spectrum.  Hence,  these  parameters, 
while  adequate  for  preliminary  analysis,  must  be  considered  as  very  rough 
characterizations . 


41.  Reflection  coefficient.  Following  notation  and  definitions  given 
in  SPM,  reflection  coefficients  x  *re  generally  defined  as  ratios  of 
reflected  to  incident  wave  heights.  Coefficients  indicating  ratios  of  energy 
would  then  vary  as  the  squares  of  height-based  coefficients  defined  here. 

Three  reflection  coefficients  are  identified.  The  first  is  a  bulk  parameter 
that  characterizes  total  reflected  to  total  Incident  energy.  These  two 
energies  are  the  volumes  under  the  reflected  and  incident  spectra,  respective¬ 
ly,  and  are  the  sources  of  the  definitions  of  and  ,  respectively. 

Because  it  is  a  wave  height  ratio  that  is  used  in  the  definition,  a  bulk 
reflection  coefficient  Xm»  i*  simply 
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The  other  two  reflection  coefficients  are  frequency  specific.  One  is  based  on 
the  square  root  of  the  ratio  of  reflected  to  incident  energy  at  the  peak 
frequency  of  the  incident  spectrum  and  suggests  how  much  of  the  energy  at  the 
frequency  spectral  peak  is  reflected.  Denoted  xp>l  >  it  is  defined  from  the 
Incident  and  reflected  frequency  spectra  given  by  Equations  10  and  11, 


\ 


Y 


7 


OF 


21 


respectively,  es 
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Because  the  peak  in  the  reflected  energy  stay  not  occur  at  the  peak  frequency 
of  the  incident  spectrum,  a  third  reflection  coefficient  is  defined  from  the 
ratio  of  frequency  spectral  energies  at  the  peak  frequency  of  the  reflected 
spectrum.  Denoted  xp,r  •  it  is  defined  by 
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42. 


adequate  for  this  preliminary  analysis. 
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Together,  these  parameters  give  a  bulk  characterization  of  some  s ?'  '  .■ 

i  r 

properties  of  a  frequency-direction  spectrum.  There  are,  of  course,  many 
other  parameters  that  can  be  defined,  but  the  present  set  is  simple  and  is  ^  f f Uh-L 

«•  Aififa3.ys,^.Restfl,ga  ^ 
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43.  A  magnetic  tape  containing  the  set  of  observed  frequency-direction 
spectra  has  been  prepared  for  copying  and  distribution.  These  spectra  are 
immediately  available  to  researchers  within  the  USACE  and  are  available  to  the 
general  public  after  June  1994.  They  can  be  obtained  by  communicating  with: 


Chief,  Field  Research  Facility 
1261  Duck  Road 
Kitty  Hawk,  NC  27949-4471 
Phone:  (919)  261-3511 
Fax:  (919)  261-4432 


44.  Appendix  A  contains  a  listing  of  the  date,  starting  time,  and  the 
characterizing  parameters  defined  previously  for  each  case.  It  is  intended  to 
be  used  as  a  kind  of  index  or  catalog  of  the  set  of  cases.  Dates  are  given 
in  the  form  yymmdd  where  yy  is  a  two-digit  year  indicator  (e.g. ,  90  means 
1990),  mm  is  the  numeric  index  of  the  calendar  month  (i.e,  01  is  January,  12 
is  December,  etc.),  and  dd  is  day  of  the  month.  All  times  are  Eastern 
Standard  Time.  A  24-hr  clock  is  used. 

45.  Some  of  the  bulk  parameters  defined  above  are  plotted  as  functions 
of  time  in  Appendix  B.  One  graph  is  shown  for  each  month  of  the  post- 
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experiment  period.  Shown  are  H*,  ,  r  depth,  #pl  , 

fora  of  Sp  r  ,  Tpl  ,  Tp<r  ,  *«,  .  XP,i  .  *nd  xp<r  ■  The  modification  to  #p >r 
involved  a  shift^of  180  deg  so  that,  in  the  figures,  this  variable  represents 
direction  toward  which  reflected  waves  propagate.  The  variable  0p>1  contin¬ 
ues  to  represent  direction  from  which  incident  waves  arrive. 


46.  The  intent  of  this  report  is  to  indicate  that  additional  frequency - 
direction  spectral  observations  have  been  added  to  the  already  extensive  FRF 
data  base  and  to  characterize  these  observations  with  some  simple  parameters. 

There  is  clearly  iuuch  more  work  to  be  done  with  this  extraordinary  data  ser, 
including  scientific  interpretation,  but  that  is  beyond  the  scope  of  this 
document.  The  descriptive  parameters  used  here  are  defined  simply  to  be 
preliminary  indicators  of  the  nature  of  the  data.  Close  examination  of  the 

u. 

parameters  displayed  in  Appendix  B  indicates  that  the  rudimentary  definitions  v  ' 

used  here  may  be  somewhat  misleading.  Hence,  there  are  a  few  caveatae  for  the 
reader . 

47.  In  plots  of  characteristic  wave  height,  it  is  clear  that  most  of 

the  total  wave  energy,  characterized  by  H^,  ,  Is  in  the  incident  spectrum 
which  is  characterized  by  .  The  so-called  reflected  energy,  indicated 

by  H»„ir  is  small,  but  is  very  seldom  zero.  It  is  possible  that  some  energy 
always  reflects  from  the  FRF  beach,  but  there  are  (at  least)  two  other  reasons 
for  this  indication  of  reflection  that  is  not  related  to  reflection  at  all. 

One  reason  is  that  there  is  usually  some  low-level,  broad  band  background 
noise  in  the  data,  which  the  MLE  algorithm  tends  to  distribute  uniformly 
around  the  horizon.  Integration  of  the  spectra  with  respect  to  direction  then 
tends  to  suggest  more  reflected  (and  incident,  but  to  a  lesser  exterrt)energy 
than  is  actually  present.  A  second  misleading  indicator  occurs  if  the  tails 
of  directional  distributions  for  incoming  waves  extend  beyond  the  arbitrarily 
assigned  90-deg  limit  for  incident  energy.  In  the  definitions  used  here,  the 
energy  represented  by  these  extended  tails  becomes  associated  with  reflection 
when  it  is  more  likely  part  of  an  onshore  distribution  of  energy  where  the 


peak  is  at  a  large,  but  less  than  90-deg,  angle  to  shore  normal.  Figure  6 
illustrates  both  these  effects.  A  low  level  of  energy  is  perceived  at  all  fgJ  - 
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directions  in  the  integrated  direction  spectrum  and  high-frequency  waves  have 
a  high  angle  of  attack,  with  energy  distribution  tails  that  exceed  90  deg. 

48.  These  effects  have  consequences  in  some  of  the  other  parameters,  as 
well.  If  there  is  significant  energy  in  distribution  tails  that  cross  or>e  of 
the  90-deg  boundaries,  the  search  algorithm  will  likely  find  one  of  those 
tails  and  the  corresponding  90-deg  boundary  as  the  peak  energy  point  for  the 
reflected  spectrum.  Hence,  there  will  be  a  number  of  cases  where  0p r  is 

90  deg  or  -90  deg.  This  is  seen  to  be  the  case  in  the  figures  in  Appendix  B. 
Though  not  as  obvious,  peak  reflected  periods  Tp  r  derived  from  peak  re¬ 
flected  frequencies  f  p  r  are  also  affected. 

49.  Finally,  reflection  coefficients,  deduced  from  the  integral  proper¬ 
ties  of  the  spectra,  are  also  somewhat  biased.  This  is  especially  noticeable 
when  low  incident  wave  conditions  exist.  When  the  low,  but  persistent,  back¬ 
ground  levels  of  the  so-called  reflected  energies  are  divided  by  low  incident 
wave  energies,  the  ratio  can  become  of  order  1,  which  is  unrealistically  high. 

50.  These  observations  do  not  entirely  preclude  the  usefulness  of  the 

parameters  defined  here.  High  energy  conditions  are  readily  identified,  as 
are  primary  directions  and  periods.  In  some  high  energy  conditions,  reflected 
wave  heights  rise  above  the  background  level  and  suggest  the  existence  of 
legitimate  reflected  energy.  The  point  here  is  that  brute  parameter  defini¬ 
tion,  even  though  reasonable  at  some  low  level  of  reliability,  can  mislead. 
Clearly,  more  elegant  definitions  are  required  for  sound  scientific  interpre¬ 
tation,  a  task  reserved  for  future  work.  - - - - — 


4V  ^ 

51.  Data  from  the  post-experiment  period  of  the  SAMSON  project  has  V 
provided  an  opportunity  to  add  to  the  FRF  data  base  of  high- resolution, 
directional -spectral  observations.  Processed  data  have  been  put  in  a  form 
that  will  be  highly  accessible  to  researchers  interested  in  nearshore  process¬ 
es.  Directional  gage,  directional  analysis  algorithms,  and  definitions  of 
characterizing  parameters  are  described.  The  parameters  used  here  are  crude, 
but  adequate  for  case  distinction  in  this  preliminary  analysis.  More  refined 
parametric  definitions  will  be  necessary  for  proper  scientific  interpretation 
of  the  data. 
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14.0  -178.0  0.093 

0.083  0.371 

0.3*7 

0.479 

13.33 

901120 

0700 

1.07 

0.97 

0.*3 

16.0  -168.0  0.103 

0.093  0.469 

0.4*1 

0.6*4 

14.12 

901120 

1000 

0.93 

0.88 

0.32 

1*.0  -172.0  0.093 

0.093  0.339 

0.396 

0.396 

13.91 

901120 

1300 

0.93 

0.91 

0.27 

18.0  -174.0  0.093 

0.083  0.300 

0.276 

0.366 

13.30 

901120 

1600 

0.99 

0.96 

0.27 

18.0  -176.0  0.093 

0.093  0.279 

0.390 

0.390 

13.32 

901120 

1900 

1.02 

0.98 

0.27 

24.0  -172.0  0.093 

0.093  0.279 

0.388 

0.388 

13.77 

901120 

2200 

1.03 

1.02 

0.26 

22.0  -176.0  0.181 

0.093  0.238 

5.200 

0.323 

13.70 

901121 

0100 

1.2* 

1.21 

0.28 

26.0  90.0  0.171 

0.093  0.228 

0.193 

0.313 

13.28 

901121 

0*00 

1.17 

1.1* 

0.26 

26.0  90.0  0.171 

0.103  0.228 

0.173 

0.289 

13.39 

901121 

0700 

1.1* 

1.11 

0.26 

28.0  -168.0  0.171 

0.103  0.233 

0.20* 

0.339 

13.95 

901121 

1000 

1.03 

1.00 

0.26 

24.0  -174.0  0.162 

0.103  0.260 

0.197 

0.3*4 

13.91 

901121 

1300 

i.a* 

1.01 

0.26 

24.0  -172.0  0.142 

0.123  0.232 

0.213 

0.322 

13.32 

901121 

1600 

1.02 

0.99 

0.2* 

20.0  -170.0  0.132 

0.132  0.239 

0.236 

0.236 

13.22 

901121 

1900 

o.oa 

0.84 

0.23 

22.0  -168.0  0.162 

0.123  0.273 

0.2*4 

0.381 

13.61 

901121 

2200 

0.77 

0.7* 

0.22 

28.0  -170.0  0.132 

0.113  0.298 

0.317 

0.378 

13.65 

901122 

0100 

0.70 

0.67 

0.20 

20.0  -164.0  0.103 

0.113  0.292 

0.308 

0.351 

13.23 

901122 

0*00 

0.67 

0.6* 

0.19 

14.0  -166.0  0.113 

0.113  0.306 

0.333 

0.353 

13.21 

901122 

0700 

0.63 

0.61 

0.22 

16.0  -164.0  0.113 

0.113  0.337 

0.333 

0.333 

13.71 

901122 

1000 

0.61 

0.38 

0.18 

16.0  -172.0  0.113 

0.113  0.315 

0.378 

0.378 

13.86 

901122 

1300 

0.36 

0.33 

0.18 

10.0  -134.0  0.113 

0.113  0.331 

0.332 

0.352 

13.36 

901122 

1600 

0.3* 

0.31 

0.16 

12.0  -172.0  0.123 

0.123  0.320 

0.35* 

0.35* 

13.11 

901122 

1900 

0.3* 

0.31 

0.18 

6.0  -150.0  0.132 

0.123  0.362 

0.3*5 

0.395 

13. *7 

901122 

2200 

0.*9 

0.*7 

0.15 

8.0  -174.0  0.132 

0.132  0.331 

0.386 

0.386 

13.68 

901123 

0100 

0.*7 

0.** 

0.16 

10.0  -132.0  0.123 

0.113  0.37* 

0.335 

0.40* 

13.32 

901123 

0*00 

o.*a 

0.*5 

0.16 

4.0  -176.0  0.132 

0.103  0.358 

0.286 

0.527 

13.11 

901123 

0700 

0.31 

0.*7 

0.19 

8.0  -152.0  0.132 

0.103  0.402 

0.338 

0.3*2 

13.31 

901123 

1000 

0 .  *7 

0.** 

0.17 

10.0  -178.0  0.113 

0.123  0.377 

0.373 

0.391 

13.80 

901123 

1300 

0 .  *3 

0.39 

0.17 

8.0  176.0  0.113 

0.113  0.437 

0.493 

0.495 

13.  *2 

901123 

1600 

0  ■  *2 

0.38 

0.18 

8.0  -174.0  0.123 

0.113  0.476 

0.407 

0.337 

13.0* 

901123 

1900 

0 .  *6 

0.*2 

0.18 

6.0  -156.0  0.123 

0.113  0.436 

0.389 

0.319 

13.27 

901123 

2200 

o.*o 

0.37 

0.13 

14.0  -176.0  0.113 

0.103  0.401 

0.33* 

0.473 

13.62 

90112* 

0100 

0.3* 

0.31 

0.13 

2.0  -174.0  0.113 

0.093  0.403 

0.38* 

0.485 

13.  *0 

90112* 

0*00 

0.29 

0.27 

0.11 

12.0  -178.0  0.093 

0.103  0 . *20 

0.40* 

0.47* 

13.06 

90112* 

0700 

0.20 

0.23 

0.12 

6.0  -176.0  0.093 

0.093  0.489 

0.39* 

0.39* 

13.3* 

90112* 

1000 

0.30 

0.27 

0.13 

4.0  -172.0  0.093 

0.093  0.481 

0.61* 

0.614 

13.75 

90112* 

1300 

0.31 

0.28 

0.1* 

20.0  -170.0  0.093 

0.093  0.302 

0.566 

0.566 

13. *9 

90112* 

1600 

0.32 

0.29 

0.13 

6.0  -172.0  0.093 

0.093  0.322 

0.529 

0.529 

13.01 

90112* 

1900 

0.3* 

0.29 

0.17 

20.0  -176.0  0.093 

0.093  0.590 

0.635 

0.633 

13.08 

90112* 

2200 

0.3* 

0.29 

0.16 

20.0  -172.0  0.093 

0.093  0.356 

0.397 

0.397 

13.47 

901123 

0100 

0.30 

0.28 

0.13 

12.0  -174.0  0.093 

0.093  0.439 

0.467 

0.467 

13.38 

901123 

0*00 

0.29 

0.26 

0.13 

16.0  -174.0  0.093 

0.093  0.303 

0.537 

0.337 

12.93 

901123 

0700 

0.23 

0.22 

0.12 

4.0  -178.0  0.093 

0.103  0.3*3 

0.369 

0.686 

12.97 

901123 

1000 

0.2* 

0.22 

0.11 

6.0  -174.0  0.103 

0.093  0.303 

0.5*3 

0.671 

13.39 

901123 

1300 

0.2* 

0.22 

0.10 

-36.0  -172.0  0.093 

0.083  0.439 

0.412 

0.626 

13.40 

901123 

1600 

0.23 

0.23 

0.10 

-38.0  -90.0  0.103 

0.103  0.426 

0.478 

0.478 

12.93 

901123 

1900 

0.26 

0.2* 

0.11 

-38.0  -90.0  0.103 

0.06*  0.467 

0.603 

0.733 

12.83 

901123 

2200 

0.28 

0.26 

0.10 

-56.0  -90.0  0.103 

0.06*  0.40* 

0.462 

0.7*1 

13.31 

901126 

0100 

0.26 

0.2* 

0.09 

-34.0  -90.0  0.162 

0.06*  0.388 

0.332 

0.710 

13.49 

901126 

0*00 

0.23 

0.23 

0.10 

-36.0  -90.0  0.132 

0.06*  Q.*3* 

0.315 

0.768 

13.11 
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Dots 

Tlaa 

1ST 

V 

^.r 

a.  .. 

Jai 

V  V 

Jfp.l 

rf.t 

Dtpth 

■ 

441126 

0700 

0.23 

0.22 

0.12 

•26.0 

-178.0 

0.064 

0.064 

0.326 

0.642 

0.642 

12.92 

901126 

1000 

0.26 

0.23 

0.13 

-34.0 

-172.0 

0.064 

0.054 

0.334 

0.729 

0.889 

11.33 

901126 

1300 

0.30 

0.26 

0.15 

-28.0 

-176.0 

0.034 

0.034 

0.363 

0.876 

0.876 

13.35 

901126 

1600 

0.32 

0.27 

0.16 

-26.0 

-174.0 

0.034 

0.034 

0.603 

0.761 

0.761 

13.13 

901126 

1900 

0.33 

0.27 

0.18 

-24.0 

-176.0 

0.034 

0.034 

0.665 

0.780 

0.780 

12.83 

901126 

2200 

0.39 

0.32 

0.23 

4.0 

-172.0 

0.034 

0.034 

0.727 

0.893 

0.895 

13.30 

901127 

0100 

0.40 

0.33 

0.22 

-6.0 

-174.0 

0.034 

0.034 

0.631 

0.742 

0.742 

13.71 

901128 

0100 

0.43 

0.40 

0.21 

6.0 

-170.0 

0.074 

0.064 

0.515 

0.684 

0.811 

13.72 

901128 

0400 

0.44 

0.40 

0.18 

2.0 

-172.0 

0.074 

0.064 

0.462 

0.344 

0.711 

13.39 

901128 

0700 

0.46 

0.42 

0.18 

14.0 

-90.0 

0.074 

0.074 

0.433 

0.382 

0.382 

13.02 

901128 

1000 

0.33 

0.51 

0.21 

-58.0 

-90.0 

0.191 

0.074 

0.404 

0.277 

0.649 

13.03 

901128 

1300 

0.35 

0.30 

0.21 

2.0 

-90.0 

0.191 

0.074 

0.409 

0.340 

0.675 

13.53 

901128 

1600 

0.38 

0.34 

0.21 

-34.0 

-90.0 

0.191 

0.064 

0.382 

0.368 

0.732 

13.48 

901128 

1900 

0.48 

0.43 

0.17 

-32.0 

-90.0 

0.074 

0.064 

0.379 

0.431 

0.717 

12.82 

901128 

2200 

0.47 

0.43 

0.19 

-36.0 

-90.0 

0.201 

0.074 

0.437 

0.341 

0.724 

12.83 

901129 

0100 

0.54 

0.51 

0.18 

-56.0 

-90.0 

0.191 

0.074 

0.353 

0.270 

0.341 

13.54 

901129 

0400 

0.56 

0.52 

0.20 

-32.0 

-90.0 

0.191 

0.074 

0.384 

0.348 

0.404 

13.69 

901129 

0700 

0.52 

0.49 

0.18 

-50.0 

-90.0 

0.171 

0.074 

0.373 

0.300 

0.586 

13.03 

901129 

1000 

1.13 

1.08 

0.34 

32.0 

90.0 

0.191 

0.240 

0.315 

0.230 

0.440 

12.80 

901129 

1300 

1.72 

1.67 

0.39 

38.0 

90.0 

0.162 

0.181 

0.233 

0.190 

0.261 

13.43 

901129 

1600 

1.68 

1.62 

0.46 

36.0 

90.0 

0.171 

0.171 

0.282 

0.263 

0.263 

13.76 

901129 

1900 

1.78 

1.71 

0.48 

34.0 

90.0 

0.171 

0.181 

0.279 

0.262 

0.297 

13.27 

901129 

2200 

1.96 

1.90 

0.50 

34.0 

90.0 

0.162 

0.171 

0.261 

0.184 

0.282 

13.03 

901130 

0100 

2.38 

2.31 

0.60 

32.0 

90.0 

0.142 

0.240 

0.261 

0.143 

0.433 

13.76 

901130 

0400 

2.13 

2.07 

0.38 

36.0 

90.0 

0.142 

0.240 

0.282 

0.171 

0.403 

14.30 

901130 

0700 

l.SS 

1.81 

0.48 

28.0 

90.0 

0.132 

0.240 

0.265 

0.197 

0.381 

13.79 

901130 

1000 

1.69 

1.63 

0.36 

30.0 

90.0 

0.142 

0.191 

0.221 

0.138 

0.273 

13.09 

901130 

1300 

1.48 

1.43 

0.31 

26.0 

90.0 

0.142 

0.240 

0.217 

0.167 

0.312 

13.37 

901130 

1600 

1.31 

1.27 

0.32 

30.0 

90.0 

0.142 

0.142 

0.236 

0.260 

0.260 

13.91 

901130 

1900 

1.13 

1.10 

0.27 

22.0 

90.0 

0.123 

0.103 

0.249 

0.190 

0.311 

13.46 

901130 

2200 

1.01 

0.99 

0.21 

22.0 

90.0 

0.123 

0.123 

0.214 

0.210 

0.210 

12.84 

901201 

0100 

0.89 

0.85 

0.25 

24.0 

-166.0 

0.142 

0.123 

0.293 

0.218 

0.339 

13.36 

901201 

0400 

0.76 

0.72 

0.26 

22.0 

-164.0 

0.142 

0.113 

0.366 

0.295 

0.342 

14.18 

901201 

0700 

0.64 

0.39 

0.23 

16.0 

-170.0 

0.132 

0.142 

0.384 

0.383 

0.425 

13.84 

901201 

1300 

0.32 

0.49 

0.17 

16.0 

180.0 

0.142 

0.103 

0.357 

0.264 

0.417 

12.88 

901201 

1900 

0.46 

0.42 

0.19 

10.0 

-176.0 

0.103 

0.103 

0.466 

0.464 

0.464 

13.44 

901202 

0100 

0.48 

0.43 

0.20 

10.0 

-178.0 

0.093 

0.093 

0.467 

0.500 

0.500 

12.88 

901202 

0700 

0.41 

0.37 

0.19 

0.0 

-174.0 

0.103 

0.093 

0.509 

0.418 

0.547 

13.97 

901202 

1300 

0.39 

0.36 

0.15 

4.0 

-176.0 

0.103 

0.103 

0.429 

0.383 

0.383 

12.66 

901202 

1900 

0.41 

0.36 

0.19 

0.0 

-170.0 

0.103 

0.093 

0.531 

0.326 

0.538 

13.71 

901203 

0100 

0.46 

0.42 

0.19 

2.0 

-176.0 

0.103 

0.083 

0.447 

0.387 

0.339 

12.88 

901203 

0700 

1.03 

1.00 

0.22 

8.0 

-90.0 

0.201 

0.083 

0.224 

0.172 

0.350 

14.33 

901203 

1300 

1.10 

1.07 

0.24 

8.0 

-90.0 

0.171 

0.093 

0.223 

0.132 

0.328 

12.84 

901203 

1900 

1.38 

1.32 

0.40 

-10.0 

-90.0 

0.132 

0.132 

0.300 

0.349 

0.349 

13.92 

901204 

0100 

1.33 

1.43 

0.49 

2.0 

-90.0 

0.132 

0.064 

0.337 

0.271 

0.334 

12.71 

901204 

0700 

1.63 

1.31 

0.39 

-2.0 

-90.0 

0.123 

0.064 

0.392 

0.372 

0.324 

14.12 

901204 

1300 

1.09 

1.02 

0.17 

0.0 

176.0 

0.074 

0.074 

0.366 

0.341 

0.341 

12.83 

901204 

1900 

0.79 

0.72 

0.32 

10.0 

-176.0 

0.074 

0.074 

0.443 

0.323 

0.523 

13.84 

901203 

0100 

0.85 

0.82 

0.23 

34.0 

90.0 

0.083 

0.083 

0.301 

0.338 

0.336 

12.93 

901203 

0700 

0.87 

0.84 

0.24 

32.0 

90.0 

0.162 

0.113 

0.293 

0.224 

0.523 

14.04 
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Data 

tlaa 

X2Z_ 

*■>.1 

“-o.r 

'5:.1  V  ‘t.- 

xP.l 

XP.‘ 

Dapcii 

901203 

1300 

0.80 

0.78 

0.20 

32.0  90.0  0.132 

0.083 

0.262 

0.21* 

0.33* 

13.03 

901203 

1900 

0.38 

0.33 

0.17 

38.0  90.0  0.171 

0.093 

0.312 

0.336 

0.396 

13.3* 

901203 

0100 

0.30 

0.  *7 

0.17 

28.0  180.0  0.083 

0.083 

0.361 

0.370 

0.370 

12.97 

901203 

0700 

o.*s 

0.*3 

0.18 

1*.0  -180.0  0.093 

0.07* 

0.389 

0.379 

0.321 

13.62 

901203 

1300 

o.*s 

0 .  *1 

0.19 

-3.0  -178.0  0.07* 

0.06* 

0.430 

0.303 

0.317 

13.1* 

901203 

1900 

0.*1 

0.37 

0.18 

-2.0  -180.0  0.074 

0.07* 

o.tas 

0.498 

0.498 

13.16 

901207 

0100 

0.38 

0.34 

0.17 

-3.0  -172.0  0.07* 

0.07* 

0.*97 

0.493 

0.493 

13.20 

901207 

0700 

0.38 

0.32 

0.17 

0.0  -180.0  0.07* 

0.07* 

0.321 

0.608 

0.608 

13.36 

901207 

1300 

0.37 

0.3* 

0.18 

•2.0  -178.0  0.083 

0.07* 

0.476 

0.366 

0.373 

13.5* 

901207 

1900 

0.*0 

0.38 

0.18 

•8.0  -172.0  0.083 

0.07* 

0.438 

0.376 

0.39* 

13.2* 

901208 

0100 

0.72 

0.71 

0.13 

22.0  90.0  0.230 

0.2*0 

0.21* 

0.139 

0.137 

13.78 

901208 

0700 

1.37 

1.3* 

0.33 

2.0  90.0  0.191 

0.191 

0.223 

0.229 

0.229 

13.39 

901208 

1300 

2.08 

2.01 

0.30 

8.0  90.0  0.171 

0.181 

0.2*9 

0.213 

0.273 

13.88 

901208 

1900 

2.03 

2.00 

0.  *2 

14.0  90.0  0.1*2 

0.132 

0.207 

0.131 

0.163 

13.32 

901209 

0100 

2.03 

2.02 

0.  *3 

22.0  90.0  0.132 

0.1*2 

0.211 

0.167 

0.191 

1*.03 

901209 

0700 

2.10 

2.03 

0.** 

22.0  90.0  0.1*2 

0.2*0 

0.217 

0.1*6 

0.368 

13.37 

901209 

1300 

1.7* 

1.89 

0.*0 

15.0  -170.0  0.123 

0.103 

0.239 

0.173 

0.31* 

13.96 

901209 

1900 

1.27 

1.23 

0.29 

18.0  -175.0  0.113 

0.093 

0.237 

0.217 

0.298 

13.12 

901210 

0100 

0.93 

0.89 

0.3* 

16.0  -172.0  0.113 

0.103 

0.382 

0.280 

0.  *90 

13.88 

901210 

0700 

0.72 

0.89 

0.22 

14.0  -178.0  0.113 

0.093 

0.318 

0.238 

0.386 

13.0* 

901210 

1300 

0.3* 

0.*9 

0.23 

13.0  -180.0  0.103 

0.093 

0.478 

0.313 

0.318 

13.37 

901210 

1900 

0.31 

0.*7 

0.20 

12.0  -90.0  0.103 

0.083 

0.426 

0.301 

0.339 

12.8* 

901211 

0100 

1.29 

1.27 

0.2* 

*0.0  90.0  0.181 

0.181 

0.193 

0.199 

0.199 

13.77 

901211 

0700 

1.23 

1.2* 

0.23 

38.0  90.0  0.132 

0.171 

0.20* 

0.1*8 

0.186 

13.19 

901211 

1300 

1.19 

1.13 

0.28 

28.0  90. 0  0.132 

0.132 

0.223 

0.213 

0.237 

13.62 

901211 

1900 

0.91 

0.89 

0.21 

24.0  -174.0  0.1*2 

0.07* 

0.238 

0.191 

0.381 

13.0* 

901212 

0100 

0.7* 

0.71 

0.21 

18.0  -170.0  0.171 

0.07* 

0.292 

0.206 

0.423 

13.72 

901212 

0700 

0.39 

0.35 

0.19 

18.0  -168.0  0.083 

0.07* 

0.3*3 

0.316 

0.337 

13.20 

901212 

1300 

0.30 

0.*8 

0.18 

10.0  -172.0  0.07* 

0.083 

0.383 

0.353 

0.49* 

13.41 

901212 

1900 

0.*3 

0.39 

0.18 

8.0  -176.0  0.07* 

0.07* 

0.409 

0.305 

0.303 

13.05 

901213 

0100 

0 .  *6 

0.*3 

0.18 

-2.0  -90.0  0.083 

0.083 

0.380 

0.468 

0.468 

13.45 

901213 

0700 

0.*1 

0.38 

0.18 

-18.0  -90.0  0.162 

0.162 

0.431 

0.433 

0.433 

13.11 

901213 

1300 

o.*o 

0.37 

0.16 

-18.0  -170.0  0.182 

0.083 

0.*3* 

0.432 

0.679 

13.11 

901213 

1900 

0.3* 

0.32 

0.13 

-28.0  -1*4.0  0.132 

0.132 

0.402 

0.391 

0.391 

13.10 

90121* 

0100 

1.11 

1.07 

0.29 

36.0  90.0  0.201 

0.220 

0.266 

0.22* 

0.26* 

13.45 

90121* 

0700 

2.17 

2.12 

0.** 

34.0  90.0  0.1*2 

0.220 

0.206 

0.130 

0.303 

13.79 

90121* 

1300 

1.38 

1.3* 

0.27 

24.0  90.0  0.1*2 

0.2*0 

0.178 

0.130 

0.231 

13.35 

90121* 

1900 

1.2* 

1.22 

0.2* 

20.0  90.0  0.123 

0.123 

0.19* 

0.18* 

0.18* 

13.  *8 

901213 

0100 

1.21 

1.19 

0.23 

16.0  90.0  0.181 

0.132 

0.193 

0.173 

0.2*9 

13.33 

901213 

0700 

1.01 

0.98 

0.23 

10.0  90.0  0.181 

0.201 

0.262 

0.2*2 

0.*13 

13.70 

901213 

1300 

1.03 

1.01 

0.22 

4.0  -90.0  0.162 

0.132 

0.219 

0.229 

0.229 

13.02 

901213 

1900 

1.09 

1.0* 

0.31 

0.0  -90.0  0.132 

0.132 

0.293 

0.279 

0.279 

13.45 

901213 

0100 

0.83 

0.79 

0.23 

-3.0  -90.0  0.1*2 

0.132 

0.288 

0.273 

0.301 

13.01 

901213 

0700 

0.83 

0.38 

0.2* 

-26.0  -123.0  0.1*2 

0.1*2 

0.409 

0.4*1 

0.4*1 

13.62 

901213 

1300 

0.38 

0.3* 

0.13 

-23.0  -124.0  0.1*2 

0.132 

0.37* 

0.313 

0.338 

12.79 

901213 

1900 

0 .  *8 

0.  *7 

0.12 

36.0  90.0  0.2*0 

0.2*0 

0.233 

0.20* 

0.20* 

13.51 

9r1217 

0100 

0.32 

0.31 

0.12 

14.0  90.0  0.2*0 

0.2*0 

0.226 

0.177 

0.177 

13.02 

(Coatlmiad) 
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(Continued) 


an  H«l«ht  Plrtetloa  _ hfllgM 


Beta 

Tlaa 

£SX_ 

ill 

X— 

*9.1 

*9.8 

Depth 

■ 

SO 1217 

0700 

0.78 

0.7* 

0.19 

3*.0 

90.0 

0.171 

0.171 

0.233 

0.237 

0.237 

13.82 

901217 

1300 

0.63 

0.62 

0.13 

2*.0 

90.0 

0.171 

0.132 

0.20* 

0.1*1 

0.233 

12.79 

901217 

1900 

0.39 

0.37 

0.16 

18.0 

90.0 

0.171 

0.152 

0.290 

0.230 

0.390 

13.48 

90121S 

0100 

0.33 

0.33 

0.1* 

10.0 

-176.0 

0.152 

0.152 

0.267 

0.2*0 

0.2*0 

12.83 

901218 

0700 

0.*6 

0.** 

0.1* 

•10.0 

-176.0 

0.132 

0.123 

0.31* 

0.275 

0.31* 

13.78 

901218 

1300 

0.*1 

0.38 

0.1* 

-12.0 

-90.0 

0.123 

0.07* 

0.333 

0.264 

0.586 

12.69 

901218 

1900 

0.33 

0.52 

0.18 

-38.0 

-90.0 

0.2*0 

0.2*0 

0.35* 

0.4*3 

0.4*3 

13.3* 

901219 

0100 

0.3* 

0.51 

0.17 

-*«.Q 

•90.0 

0.123 

0.093 

0.326 

0.223 

0.363 

12.68 

901219 

0700 

0.33 

0.30 

0.19 

-38.0 

-90.0 

0.103 

0.103 

0.375 

0.473 

0.473 

13.73 

901219 

1300 

0.31 

0.*9 

0.13 

-10.0 

•90 . 0 

0.103 

0.093 

0.26* 

0.210 

0.318 

12.90 

901219 

1900 

1.03 

1.02 

0.23 

8.0 

90.0 

0.201 

0.181 

0.222 

0.135 

0.239 

13.68 

901220 

0100 

1.17 

1.13 

0.21 

16.0 

90.0 

0.181 

0.181 

0.186 

0.183 

0.183 

13.07 

901220 

0700 

l.*6 

1.** 

0.28 

16.0 

90.0 

0.171 

0.181 

0.197 

0.165 

0.189 

13.99 

901220 

1300 

1.53 

1.31 

0.2* 

12.0 

90.0 

0.162 

0.181 

0.162 

0.130 

0.183 

13.16 

901220 

1900 

1.38 

1.53 

0.32 

8.0 

90.0 

0.132 

0.1*2 

0.203 

0.158 

0.2** 

13.71 

901221 

0100 

1.36 

1.53 

0.30 

12.0 

-90.0 

0.152 

0.123 

0.196 

0.131 

0.198 

13.17 

901221 

0700 

1.32 

l.*6 

0.35 

0.0 

-17*. 0 

0.1*2 

0.132 

0.237 

0.213 

0.233 

13.83 

901221 

1300 

1.33 

1.52 

0.31 

-6.0 

-90.0 

0.123 

0.123 

0.203 

0.130 

0.130 

13.16 

901222 

0100 

1.33 

1.32 

0.30 

-14.0 

-90.0 

0.123 

0.113 

0.227 

0.139 

0.189 

13.22 

901222 

0700 

1.3* 

1.30 

0.31 

-*.0 

-90.0 

0.123 

0.103 

0.239 

0.191 

0.293 

13.38 

901222 

1300 

1.05 

1.01 

0.27 

-8.0 

-90.0 

0.123 

0.1*2 

0.263 

0.206 

0.270 

13.21 

901222 

1900 

1.03 

1.00 

0.23 

-*.0 

-90.0 

0.1*2 

0.123 

0.231 

0.223 

0.2*2 

13.3* 

901223 

0100 

0.89 

0.86 

0.26 

-2.0 

-90.0 

0.113 

0.083 

0.306 

0.23* 

0.439 

13.28 

901223 

0700 

1.11 

1.06 

0.31 

0.0 

-90.0 

0.132 

0.132 

0.29* 

0.272 

0.272 

13.2* 

901223 

1300 

1.01 

0.96 

0.31 

0.0 

-90.0 

0.1*2 

0.152 

0.320 

0.308 

0.330 

13.18 

901223 

1900 

0.91 

0.88 

0.2* 

-6.0 

-90.0 

0.093 

0.093 

0.270 

0.2*3 

0.2*3 

13.00 

90122* 

0100 

0.82 

0.76 

0.29 

-2.0 

•1*0.0 

0.093 

0.093 

0.383 

0.328 

0.328 

13.27 

90122* 

0700 

0.91 

0.86 

0.29 

-2*.0 

-90.0 

0.093 

0.093 

0.332 

0.283 

0.283 

12.93 

90122* 

1300 

0.93 

0.88 

0.28 

-1*.0 

90.0 

0.093 

0.093 

0.322 

0.238 

0.238 

13.30 

90122* 

1900 

1.87 

1.8* 

0.3* 

30.0 

90.0 

0.152 

0.210 

0.187 

0.106 

0.273 

13.18 

901223 

0100 

1.56 

1.32 

0.35 

30.0 

90.0 

0.132 

0.2*0 

0.231 

0.131 

0.350 

13.73 

901223 

0700 

1.18 

1.15 

0.23 

30.0 

90.0 

0.181 

0.2*0 

0.221 

0.186 

0.290 

13.07 

901223 

1300 

1.09 

1.06 

0.27 

32.0 

90.0 

0.093 

0.093 

0.233 

0.289 

0.289 

13.49 

901223 

1900 

1.03 

1.02 

0.23 

•8.0 

90.0 

0.083 

0.083 

0.2*1 

0.186 

0.186 

12.89 

901228 

0100 

1.02 

0.96 

0.3* 

-10.3 

-168.0 

0.083 

0.083 

0.338 

0.390 

0.390 

13.77 

901228 

0700 

0.9* 

0.90 

0.27 

-1*.0 

-172.0 

0.083 

0.083 

0.293 

0.303 

0.303 

13.07 

901226 

1300 

1.17 

1.1* 

0.29 

6.0 

90.0 

0.083 

0.083 

0.233 

0.293 

0.293 

13.69 

901226 

1900 

1.21 

1.18 

0.26 

6.0 

90.0 

0.191 

0.201 

0.218 

0.191 

0.207 

12.97 

901227 

0100 

1.38 

1.33 

0.30 

16.0 

90.0 

0.191 

0.093 

0.220 

0.175 

0.289 

13.97 

901227 

0700 

1.67 

1.63 

0.3* 

16.0 

90.0 

0.181 

0.191 

0.209 

0.17* 

0.211 

13.19 

901227 

1300 

1.91 

1.87 

0.38 

12.0 

90.0 

0.162 

0.162 

0.202 

0.170 

0.170 

13.7* 

901227 

1900 

1.79 

1.77 

0.29 

12.0 

90.0 

0.152 

0.152 

0.166 

0.1*3 

0.143 

13.06 

901228 

0100 

1.90 

1.84 

0.** 

1*.0 

-90.0 

0.132 

0.123 

0.237 

0.178 

0.209 

13.97 

901228 

0700 

1.60 

1.36 

0.39 

12.0 

-172.0 

0.113 

0.123 

0.2*8 

0.212 

0.223 

13.33 

901228 

1300 

1.29 

1.23 

0.39 

1*.0 

-170.0 

0.113 

0.103 

0.317 

0.270 

0.300 

13.67 

901228 

1900 

1.01 

0.98 

0.27 

14.0 

180.0 

0.113 

0.103 

0.279 

0.2*3 

0.321 

13.17 

901229 

0100 

1.0* 

0.97 

0.37 

10.0 

-180.0 

0.103 

0.103 

0.379 

0.432 

0.452 

13.82 
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(Conrlmud) 


an  Dir«etion_  ..famtag  _ gtilasUra. 


Oaca 

Tiaa 

1ST 

■■o.i 

tmo.t 

■ 

*mo 

*P.i 

*P.r 

Dapth 

_a_ 

901229 

0700 

1.08 

1.02 

0.33 

16.0 

180.0 

0.103 

0.103 

0.361 

0.317 

0.317 

13.60 

901229 

1300 

1.16 

1.11 

0.32 

10.0 

-176.0 

0.093 

0.093 

0.287 

0.236 

0.236 

13.63 

901229 

1900 

1.06 

1.01 

0.32 

2.0 

-176.0 

0.093 

0.093 

0.316 

0.311 

0.311 

13.31 

901230 

0100 

1.00 

0.93 

0.31 

8.0 

-90.0 

0.093 

0.083 

0.323 

0.232 

0.600 

13.38 

901230 

0700 

0.89 

0.81 

0.37 

*8.0 

-90.0 

0.093 

0.093 

0.637 

0.323 

0.323 

13.66 

901230 

1300 

0.83 

0.78 

0.27 

-2.0 

•90.0 

0.093 

0.093 

0.366 

0.326 

0.326 

12.83 

901230 

1900 

0.73 

0.69 

0.29 

•10.0 

-132.0 

0.093 

0.093 

0.626 

0.373 

0.373 

13.23 

901231 

0100 

0.67 

0.62 

0.23 

-12.0 

176.0 

0.093 

0.093 

0.376 

0.362 

0.362 

12.72 

901231 

1300 

1.67 

1.66 

0.33 

36.0 

90.0 

0.162 

0.260 

0.199 

0.127 

0.287 

12.82 

901231 

1900 

1.66 

1.62 

0.33 

22.0 

90.0 

0.132 

0.210 

0.230 

0.131 

0.277 

13.79 

910101 

0100 

1.38 

1.33 

0.31 

20.0 

90.0 

0.201 

0.260 

0.226 

0.227 

0.317 

12.91 

910101 

0700 

1.61 

1.37 

0.33 

16.0 

90.0 

0.181 

0.171 

0.233 

0.199 

0.276 

16.26 

910101 

1300 

1.16 

1.13 

0.23 

12.0 

90.0 

0.162 

0.093 

0.207 

0.130 

0.229 

12.87 

910101 

1900 

1.18 

1.16 

0.29 

10.0 

90.0 

0.181 

0.093 

0.236 

0.202 

0.331 

13.93 

910102 

0100 

1.03 

1.00 

0.26 

6.0 

-90.0 

0.093 

0.093 

0.233 

0.231 

9.231 

12.82 

910102 

0700 

1.13 

1.10 

0.36 

6.0 

-166.0 

0.162 

0.162 

0.303 

0.331 

0.331 

16.28 

910102 

1300 

1.16 

1.13 

0.28 

18.0 

90.0 

0.230 

0.230 

0.232 

0.281 

0.281 

12.89 

910102 

1900 

1.27 

1.22 

0.36 

16.0 

-176.0 

0.171 

0.181 

0.293 

0.268 

0.293 

13.96 

910103 

0100 

0.93 

0.91 

0.20 

16.0 

-90.0 

0.171 

0.181 

0.223 

0.187 

0.263 

12.83 

910103 

0700 

0.78 

0.76 

0.23 

16.0 

-176.0 

0.171 

0.103 

0.361 

0.286 

0.660 

16.08 

910103 

1300 

0.63 

0.61 

0.18 

6.0 

-90.0 

0.123 

0.103 

0.288 

0.238 

0.303 

12.92 

910103 

1900 

0.73 

0.69 

0.22 

8.0 

-178.0 

0.123 

0.103 

0.311 

0.303 

0.380 

13.81 

910109 

0100 

0.37 

0.36 

0.18 

6.0 

-180.0 

0.103 

0.103 

0.326 

0.327 

0.327 

13.07 

910109 

0700 

1.22 

1.19 

0.26 

36.0 

90.0 

0.201 

0.191 

0.220 

0.163 

0.213 

13.99 

910109 

1300 

1.36 

1.32 

0.27 

36.0 

90.0 

0.171 

0.171 

0.203 

0.176 

0.176 

13.26 

910109 

1900 

1.03 

1.03 

0.22 

20.0 

90.0 

0.181 

0.181 

0.209 

0.183 

0.183 

13.67 

910103 

0100 

1.03 

1.03 

0.20 

20.0 

90.0 

0.181 

0.260 

0.197 

0.136 

0.196 

13.60 

910103 

0700 

1.18 

1.13 

0.26 

16.0 

90.0 

0.201 

0.201 

0.209 

0.199 

0.199 

13.63 

910103 

1300 

1.03 

1.00 

0.26 

8.0 

-90.0 

0.181 

0.181 

0.236 

0.213 

0.213 

13.60 

910103 

1900 

0.97 

0.96 

0.26 

2.0 

-90.0 

0.132 

0.132 

0.231 

0.261 

0.261 

13.39 

910106 

0100 

0.73 

0.71 

0.26 

6.0 

-176.0 

0.132 

0.132 

0.360 

0.367 

0.367 

13.36 

910106 

0700 

0.61 

0.38 

0.19 

8.0 

-90.0 

0.162 

0.162 

0.332 

0.310 

0.310 

13.23 

910106 

1300 

0.67 

0.66 

0.18 

-8.0 

-178.0 

0.103 

0.113 

0.613 

0.608 

0.638 

13.66 

910106 

1900 

0.30 

0.66 

0.18 

-2.0 

-90.0 

0.113 

0.113 

0.396 

0.637 

0.637 

13.10 

910107 

0100 

0.32 

0.68 

0.21 

-8.0 

-176.0 

0.113 

0.113 

0.627 

0.636 

0.636 

13.70 

910107 

0700 

0.30 

0.66 

0.19 

0.0 

-178.0 

0.113 

0.103 

0.613 

0.606 

0.663 

13.23 

910107 

1300 

1.98 

1.93 

0.66 

36.0 

90.0 

0.132 

0.162 

0.228 

0.170 

0.196 

13.87 

910107 

1900 

2.92 

2.87 

0.32 

10.0 

90.0 

0.123 

0.260 

0.179 

0.130 

0.336 

13.67 

m 

o 

H 

O 

H 

« 

0100 

3.00 

2.93 

0.63 

12.0 

90.0 

0.123 

0.113 

0.213 

0.177 

0.198 

16.22 

910108 

0700 

2.73 

2.70 

0.69 

10.0 

90.0 

0.113 

0.260 

0.183 

0.136 

0.327 

13.30 

910108 

1300 

3.01 

2.97 

0.36 

16.0 

90.0 

0.123 

0.132 

0.182 

0.163 

0.169 

16.00 

910108 

1900 

2.86 

2.79 

0.68 

6.0 

90.0 

0.093 

0.260 

0.172 

0.131 

0.396 

13.39 

910109 

0100 

3.11 

3.03 

0.39 

0.0 

90.0 

0.093 

0.093 

0.193 

0.137 

0.137 

16.19 

910109 

0700 

2.69 

2.66 

0.36 

2.0 

90.0 

0.083 

0.260 

0.207 

0.132 

0.309 

13.68 

910109 

1300 

2.61 

2.33 

0.32 

0.0 

90.0 

0.083 

0.083 

0.203 

0.126 

0.126 

13.91 

910109 

1900 

1.83 

1.78 

0.66 

-16.0 

-180.0 

0.083 

0.083 

0.266 

0.213 

0.213 

13.33 

910110 

0100 

1.69 

1.60 

0.32 

22.0 

-176.0 

0.093 

0.093 

0.327 

0.333 

0.333 

16.06 
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W«'r»  mm  Direction  Freau.ncT  tofl.ctlm 


Oat* 

Tlaa 

US- 

w 

^o.r 

_!Sii  Jlil  Jlii  Jiil 

*9.1 

Depth 

■ 

*10110 

0700 

1.31 

1.44 

0.46 

18.0  -176.0  0.103 

0.103 

0.321 

0.368 

0.368 

13.30 

910110 

1300 

1.49 

1.43 

0.44 

22.0  -172.0  0.093 

0.093 

0.309 

0.482 

0.482 

13.77 

910110 

1900 

1.13 

1.09 

0.37 

16.0  -174.0  0.093 

0.0*3 

0.337 

0.433 

0.433 

13.32 

910111 

0100 

1.47 

1.42 

0.39 

20.0  90.0  0.191 

0.103 

0.276 

0.264 

0.418 

13.91 

910111 

0700 

1.64 

1.39 

0.41 

16.0  -180.0  0.132 

0.152 

0.261 

0.243 

0.243 

13.33 

910111 

1300 

1.36 

1.30 

0.43 

12.0  178.0  0.142 

0.132 

0.302 

0.234 

0.304 

13.33 

910111 

1900 

2.00 

1.92 

0.36 

0.0  -90.0  0.123 

0.123 

0.293 

0.282 

0.282 

13.31 

910112 

0100 

2.39 

2.28 

0.70 

0.0  -90.0  0.113 

0.113 

0.303 

0.313 

0.313 

13.57 

910112 

0700 

1.66 

1.32 

0.66 

0.0  -172.0  0.103 

0.103 

0.431 

0.490 

0.490 

13.37 

910112 

1300 

1.41 

1.27 

0.60 

0.0  -174.0  0.113 

0.103 

0.469 

0.491 

0.349 

13.37 

910112 

1900 

1.16 

1.04 

0.32 

0.0  -172.0  0.113 

0.113 

0.501 

0.494 

0.494 

13.44 

910113 

0100 

1.16 

1.03 

0.49 

14.0  -172.0  0.113 

0.113 

0.463 

0.531 

0.531 

13.54 

910113 

0700 

1.43 

1.33 

0.46 

26.0  90.0  0.181 

0.181 

0.342 

0.313 

0.313 

13.83 

910113 

1300 

1.14 

1.09 

0.32 

32.0  90.0  0.171 

0.181 

0.289 

0.254 

0.321 

13.32 

910113 

1900 

1.08 

1.02 

0.33 

30.0  166.0  0.171 

0.171 

0.342 

0.362 

0.362 

13.61 

910114 

0100 

0.87 

0.83 

0.23 

28.0  -174.0  0.171 

0.171 

0.300 

0.283 

0.283 

13.41 

910114 

0700 

0.74 

0.69 

0.27 

30.0  166.0  0.171 

0.171 

0.398 

0.468 

0.468 

13.91 

910114 

1300 

0.33 

0.31 

0.20 

-2.0  -178.0  0.103 

0.064 

0.393 

0.368 

0.649 

13.07 

910114 

1900 

0.33 

0.48 

0.22 

0.0  -174.0  0.103 

0.064 

0.485 

0.390 

0.733 

13.49 

910113 

0100 

0.49 

0.43 

0.21 

-4.0  -174.0  0.064 

0.064 

0.473 

0.393 

0.593 

12  99 

910113 

0700 

0.43 

0.41 

0.19 

•8.0  -170.0  0.064 

0.064 

0.479 

0.487 

0.487 

13.78 

910113 

1300 

0.46 

0.41 

0.21 

-4.0  -176.0  0.064 

0.064 

0.524 

0.708 

0.708 

12.89 

910113 

1900 

0.46 

0.41 

0.20 

-12.0  -174.0  0.123 

0.064 

0.492 

0.308 

0.720 

13.63 

910116 

0100 

0.43 

0.41 

0.18 

•8.0  -174.0  0.074 

0.074 

0.433 

0.530 

0.550 

13.01 

910116 

0700 

1.73 

1.62 

0.63 

-10.0  -126.0  0.123 

0.113 

0.401 

0.346 

0.428 

13.90 

910116 

1300 

1.42 

1.33 

0.49 

-8.0  -110.0  0.113 

0.103 

0.369 

0.330 

0.364 

12.86 

910116 

1900 

1.08 

0.97 

0.47 

•12.0  -128.0  0.123 

0.123 

0.484 

0.300 

0.500 

13.73 

910117 

0100 

0.78 

0.72 

0.30 

-16.0  -112.0  0.113 

0.113 

0.412 

0.415 

0.415 

12.94 

910117 

0700 

0.63 

0.36 

0.33 

-14.0  -142.0  0.103 

0.103 

0.586 

0.653 

0.653 

13.94 

910117 

1300 

0.43 

0.41 

0.18 

-12.0  -178.0  0.113 

0.113 

0.450 

0.431 

0.431 

12.88 

910117 

1900 

0.43 

0.39 

0.17 

-6.0  -174.0  0.113 

0.113 

0.447 

0.484 

0.484 

13.74 

910118 

0100 

0.36 

0.33 

0.13 

-2.0  -180.0  0.113 

0.103 

0.446 

0.405 

0.310 

12.86 

910118 

0700 

0.37 

0.34 

0.13 

32.0  -180.0  0.103 

0.103 

0.387 

0.422 

0.422 

13.78 

910118 

1300 

0.30 

0.27 

0.12 

32.0  178.0  0.093 

0.093 

0.436 

0.483 

0.483 

12.74 

910118 

1900 

0.23 

0.22 

0.11 

•2.0  180.0  0.103 

0.034 

0.301 

0.561 

0.861 

13.37 

910119 

0100 

0.18 

0.17 

0.08 

•10.0  -178.0  0.103 

0.054 

0.498 

0.456 

0.698 

12.76 

910119 

0700 

0.21 

0.18 

0.10 

-26.0  -176.0  0.162 

0.034 

0.553 

0.371 

0.909 

13.63 

910119 

1300 

0.24 

0.23 

0.09 

14.0  180.0  0.103 

0.054 

0.408 

0.433 

0.692 

12.78 

910119 

1900 

0.28 

0.26 

0.11 

14.0  -180.0  0.162 

0.103 

0.414 

0.233 

0.530 

13.41 

910120 

0100 

0.23 

0.23 

0.10 

0.0  -180.0  0.103 

0.054 

0.447 

0.439 

0.908 

12.87 

910120 

0700 

0.46 

0.44 

0.13 

32.0  90.0  0.240 

0.240 

0.355 

0.251 

0.251 

13.57 

910120 

1300 

0.78 

0.73 

0.19 

20.0  90.0  0.220 

0.220 

0.260 

0.206 

0.206 

13.00 

910120 

1900 

0.63 

0.38 

0.24 

20.0  -136.0  0.220 

0.083 

0.416 

0.221 

0.638 

13.49 

910121 

0100 

0.34 

0.47 

0.26 

-4.0  -138.0  0.093 

0.093 

0.562 

0.397 

0.597 

13.25 

910121 

0700 

0.46 

0.41 

0.22 

•6.0  -140.0  0.093 

0.103 

0.548 

0.605 

0.624 

13.42 

910121 

1300 

0.33 

0.49 

0.21 

60.0  90.0  0.230 

0.230 

0.430 

0.403 

0.403 

13.30 

910121 

1900 

1.43 

1.37 

0.48 

38.0  90.0  0.162 

0.171 

0.331 

0.300 

0.379 

13.68 
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(Coesisuad) 


_ Hast  a«itht  Dlr«ctlan  fnmmT  _ K«fl«ctlOft 

^  ^o.l  *bo.«  *|,1  *p.r  *1.1  *g^r  *m>  *p.l  *p>r  O^Oj 


9X0122 

0100 

1.35 

1.27 

0.46 

36.0  90.0  0.162 

0.171 

0.366 

0.286 

0.402 

13.84 

910122 

0700 

1.59 

1.50 

0.53 

30.0  90.0  0.181 

0.201 

0.336 

0.294 

0.338 

13.67 

910122 

1300 

1.65 

1.55 

0.58 

32.0  90.0  0.152 

0.142 

0.373 

0.294 

0.434 

13.72 

910122 

1900 

1.30 

1.21 

0.47 

26.0  176.0  0.123 

0.132 

0.389 

0.444 

0.467 

13.47 

910123 

0100 

0.94 

0.88 

0.34 

32.0  164.0  0.142 

0.113 

0.381 

0.388 

0.532 

13.90 

910123 

0700 

0.70 

0.65 

0.27 

22.0  176.0  0.113 

0.113 

0.419 

0.490 

0.490 

13.31 

910123 

1300 

0.55 

0.50 

0.25 

14.0  178.0  0.113 

0.11.) 

0.494 

0.584 

0.564 

13.39 

910123 

1900 

0.52 

0.48 

0.21 

10.0  178.0  0.113 

0.103 

0.432 

0.384 

0.461 

13.00 

910124 

0100 

0.39 

0.34 

0.19 

8.0  -172.0  0.103 

0.113 

0.349 

0.509 

0.64* 

13.71 

910124 

0700 

0.38 

0.32 

0.16 

2.0  -176.0  0.093 

0.093 

0.511 

0.470 

0.470 

12.92 

910124 

1300 

0.36 

0.32 

0.16 

0.0  -174.0  0.093 

0.093 

0.306 

0.476 

0.476 

13.45 

910124 

1900 

0.36 

0.32 

0.16 

-4.0  -176.0  0.093 

0.093 

0.504 

0.492 

0.492 

12.79 

910125 

0100 

0.36 

0.32 

0.15 

-2.0  -172.0  0.093 

0.093 

0.472 

0.300 

0.500 

13.86 

910125 

0700 

1.84 

1.79 

0.44 

40.0  90.0  0.162 

0.181 

0.243 

0.154 

0.269 

13.22 

910125 

1300 

1.57 

1.49 

0.49 

32.0  90.0  0.181 

0.181 

0.328 

0.370 

0.370 

13.74 

910125 

1900 

1.10 

1.05 

0.33 

22.0  90.0  0.132 

0.132 

0.312 

0.324 

0.324 

12.96 

910126 

0100 

0.90 

0.85 

0.31 

22.0  170.0  0.132 

0.142 

0.362 

0.373 

0.419 

13.91 

910126 

0700 

0.75 

0.72 

0.22 

10.0  180.0  0.123 

0.113 

0.313 

0.309 

0.363 

13.18 

910126 

1300 

0.68 

0.64 

0.23 

10.0  -180.0  0.152 

0.123 

0.359 

0.381 

0.466 

13.33 

910126 

1900 

0.55 

0.51 

0.19 

2.0  -180.0  0.132 

0.132 

0.373 

0.377 

0.377 

12.93 

910127 

0100 

0.45 

0.41 

0.17 

0.0  -176.0  0.132 

0.142 

0.399 

0.435 

0.467 

13.66 

910127 

0700 

0.34 

0.31 

0.14 

-8.0  -174.0  0.132 

0.123 

0.466 

0.404 

0.448 

13.27 

910127 

1300 

0.33 

0.30 

0.14 

-6.0  -178.0  0.142 

0.142 

0.463 

0.487 

0.487 

13.30 

910127 

1900 

0.36 

0.33 

0.14 

-16.0  -176.0  0.132 

0.132 

0.434 

0.435 

0.455  13.11 

910128 

0100 

0.49 

0.46 

0.18 

-14.0  -90.0  0.132 

0.142 

0.385 

0.406 

0.427 

13.35 

910128 

0700 

0.63 

0.57 

0.27 

-6.0  -90.0  0.132 

0.132 

0.462 

0.346 

0.546 

13.53 

910128 

1300 

0.51 

0.48 

0.18 

-12.0  -90.0  0.132 

0.132 

0.376 

0.413 

0.413 

13.02 

910128 

1900 

0.64 

0.55 

0.32 

-12.0  -130.0  0.113 

0.113 

0.580 

0.621 

0.621 

13.38 

910129 

0100 

0.62 

0.57 

0.25 

-8.0  -108.0  0.123 

0.113 

0.448 

0.416 

0.503 

13.14 

910129 

0700 

0.61 

0.54 

0.28 

-12.0  -132.0  0.113 

0.113 

0.316 

0.537 

0.537 

13.95 

910129 

1300 

0.66 

Q.65 

0.21 

-10.0  -90.0  0.123 

0.113 

0.321 

0.303 

0.385 

12.94 

910129 

1900 

0.69 

0.65 

0.25 

-4.0  -122.0  0.132 

0.113 

0.381 

0.421 

0.500 

13.73 

910130 

0100 

0.84 

0.80 

0.26 

-6.0  -90.0  0.162 

0.171 

0.329 

0.331 

0.379 

12.93 

910130 

0700 

0.71 

0.65 

0.28 

-6.0  -172.0  0.132 

0.152 

o.4n 

0.447 

0.447 

14.08 

910201 

1300 

0.81 

0.80 

0.15 

26.0  90.0  0.152 

0.152 

0.185 

0.171 

0.171 

12.66 

910201 

1900 

0.66 

0.63 

0.19 

24.0  166.0  0.162 

0.113 

0.299 

0.269 

0.331 

13.70 

910202 

0100 

0.54 

0.53 

0.13 

6.0  90.0  0.152 

0.113 

0.255 

0.213 

0.313 

12.75 

910202 

0700 

0.55 

0.52 

0.19 

2.0  176.0  0.103 

0.103 

0.359 

0.439 

0.439 

13.71 

910202 

1300 

0.52 

0.49 

0.17 

-2.0  -104.0  0.103 

0.103 

0.335 

0.327 

0.327 

12.78 

910202 

1900 

0.54 

0.50 

0.21 

-12.0  -110.0  0.103 

0.113 

0.418 

0.481 

0.502 

13.54 

910203 

0100 

0.50 

0.46 

0.20 

•12.0  -112.0  0.103 

0.103 

0.444 

0.468 

0.468 

12.99 

910203 

0700 

0.51 

0.46 

0.22 

-6.0  -108.0  0.103 

0.103 

0.470 

0.521 

0.521 

13.44 

910203 

1300 

0.46 

0.41 

0.20 

-10.0  -116.0  0.103 

0.093 

0.476 

0.433 

0.328 

12.92 

910203 

1900 

0.54 

0.48 

0.24 

-16.0  -112.0  0.113 

0.113 

0.498 

0.513 

0.513 

13.28 

910204 

0100 

0.46 

0.40 

0.23 

-20.0  -136.0  0.103 

0.103 

0.381 

0.632 

0.632 

13.17 

910204 

0700 

0.45 

0.41 

0.20 

-10.0  -110.0  0.113 

0.103 

0.490 

0.469 

0.520 

13.22 

910204 

1300 

0.37 

0.32 

0.17 

-10.0  -136.0  0.103 

0.103 

0.327 

0.532 

0.332  13.12 

(Conelauad) 
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•1020*  1900  0.41  0.37  0.1*  -16.0  -100.0  0.113  0.113  0.*tl  0.426  0.42*  13.13 

910203  0100  0.41  0.33  0.20  -24.0  -134.0  0-113  0.103  0.337  0.331  0.041  13.47 
•10203  0700  0.43  0.3S  0.19  -21.0  -110.0  0.103  0.113  0.30*  0.471  0.333  13.21 
910203  1300  0.40  0.42  0.23  -14.0  -130.0  0.113  0.123  0.337  0.407  0.010  13.44 
910203  1900  0.31  0.40  0.23  -14.0  -132.0  0.113  0.113  0.499  0.334  0.334  13.20 

910200  0100  0.00  0.32  0.30  -10.0  -130.0  0.113  0.113  0.373  0.339  0.339  13.74 
910200  0700  0.03  0.30  0.29  -20.0  -132.0  0.113  0.113  0.314  0.309  0.309  13.10 
910200  1300  0.30  0.49  0.20  -10.0  -130.0  0.113  0.113  0.309  0.391  0.391  13.31 
910200  1900  0.33  0.40  0.24  -14.0  -120.0  0.113  0.123  0.311  0.430  0.337  13.12 

910207  0100  0.31  0.43  0.23  -10.0  -130.0  0.113  0.113  0.309  0.033  0.033  13.77 
910207  0700  0.31  0.43  0.24  -12.0  -134.0  0.123  0.123  0.322  0.303  0.303  13. IB 
910207  1300  0.34  0.47  0.23  -12.0  -140.0  0.113  0.113  0.334  0.3B3  0.3S3  13.32 
910207  1900  0.30  0.30  0.23  -0.0  -13B.0  0.123  0.113  0.304  0.4B0  0.341  13.17 

910208  0100  0.60  0.60  0.27  4.0  -144.0  0.113  0.113  0.446  0.312  0.312  13.92 
910208  0700  1.99  1.90  0.39  30.0  90.0  0.171  0.171  0.308  0.330  0.330  13.42 

910209  1300  1.33  1.26  0.43  18.0  90.0  0.210  0.123  0.340  0.336  0.433  13.33 
910209  1900  0.90  0.83  0.29  16.0  180.0  0.123  0.142  0.3*1  0.342  0.426  13.22 

910210  0100  0.71  0.66  0.28  10.0  178.0  0.123  0.123  0.386  0.437  0.437  13.67 
910210  0700  0.67  0.63  0.24  2.0  -178.0  0.113  0.113  0.384  0.389  0.389  13.41 

910211  1300  0.64  0.63  0.13  14.0  90.0  0.240  0.230  0.244  0.204  0.206  13.34 
910211  1900  0.36  0.34  0.16  8.0  -180.0  0.162  0.162  0.294  0.280  0.280  13.48 

910212  0100  1.11  1.07  0.31  26.0  90.0  0.210  0.210  0.286  0.277  0.277  13.48 
910212  0700  1.69  1.64  0.42  30.0  90.0  0.132  0.171  0.238  0.212  0.281  13.84 
910212  1300  1.16  1.13  0.27  28.0  90.0  0.162  0.171  0.244  0.219  0.270  13.29 
910212  1900  0.91  0.87  0.26  26.0  170.0  0.132  0.1*2  0.292  0.216  0.309  13.38 

910213  0100  0.64  0.61  0.18  16.0  174.0  0.123  0.123  0.302  0.322  0.322  13.13 
910213  0700  0.49  0.46  0.16  0.0  174.0  0.132  0.103  0.344  0.277  0.330  13.39 
910213  1300  0.33  0.33  0.11  -18.0  -90.0  0.201  0.103  0.349  0.320  0.336  12.81 
910213  1900  0.27  0.23  0.09  -20.0  -90.0  0.103  0.0S3  0.367  0.361  0.319  13.49 

910214  0100  0.41  0.38  0.13  -20.0  -106.0  0.132  0.132  0.399  0.408  0.408  12.93 
910214  0700  0.62  0.33  0.30  -30.0  -124.0  0.132  0.093  0.333  0.324  0.683  13.74 
910214  1300  0.66  0.61  0.27  -22.0  -90.0  0.083  0.083  0.431  0.467  0.467  12.81 
910214  1900  0.63  0.36  0.29  -20.0  -132.0  0.093  0.103  0.327  0.339  0.608  13.77 

910213  0100  0.33  0.49  0.22  -12.0  -104.0  0.103  0.093  0.432  0.419  0.305  12.94 
910213  0700  0.41  0.36  0.19  -8.0  -152.0  0.083  0.083  0.513  0.303  0.303  13.88 
910213  1300  0.39  0.36  0.13  -16.0  -102.0  0.093  0.093  0.420  0.49*  0.49*  12.81 
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910217  0700  0.29  0.26  0.13  2.0  168.0  0.083  0.083  0.496  0.350  0.330  13.60 
910217  1300  0.23  0.23  0.09  -12.0  176.0  0.093  0.064  0.392  0.303  0.581  12.52 
910217  1900  0.29  0.26  0.13  -30.0  168.0  0.093  0.083  0.483  0.394  0.731  13.37 
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0.083 

0.083 

0.391 

0.341 

0.341 

12.98 

910319 

1900 

0.84 

0.76 

0.37 

2.0 

90.0 

0.093 

0.093 

0.497 

0.311 

0.311 

14.04 

910320 

0100 

0.62 

0.36 

0.23 

10.0 

90.0 

0.083 

0.083 

0.433 

0.463 

0.463 

13.16 

910320 

0700 

0.31 

0.46 

0.23 

8.0 

174.0 

0.083 

0.093 

0 .300 

0.301 

0.623 

13.37 

910320 

1300 

0.44 

0.40 

0.18 

6.0 

138.0 

0.083 

0.083 

0.447 

0.439 

0.439 

13.02 

910320 

1900 

0.46 

0.41 

0.20 

12.0 

170.0 

0.083 

0.083 

0.486 

0.493 

0.493 

13.66 

910321 

0100 

0.48 

0.43 

0.20 

8.0 

176.0 

0.083 

0.074 

0.472 

0.403 

0.463 

13.34 

910321 

0700 

0.32 

0.46 

0.24 

8.0 

174.0 

0.074 

0.074 

0.316 

0.604 

0.604 

13.33 

910321 

1300 

0.31 

0.46 

0.22 

4.0 

-180.0 

0.074 

0.074 

0.467 

0.469 

0.469 

13-21 

910321 

1900 

0.31 

0.46 

0.23 

2.0 

162.0 

0.083 

0.083 

0.492 

0.481 

0.481 

13.33 

910322 

0100 

0.48 

0.43 

0.23 

4.0 

-166.0 

0.083 

0.074 

0.344 

0.406 

0.606 

13.52 

910322 

0700 

0.43 

0.41 

0.18 

0.0 

136.0 

0.074 

0.083 

0.446 

0.398 

0.430 

13.11 

910323 

1300 

0.82 

0.76 

0.31 

18.0 

132.0 

0.132 

0.132 

0.410 

0.393 

0.393 

13.74 

910323 

1900 

0.99 

0.93 

0.31 

16.0 

-90.0 

0.113 

0.113 

0.333 

0.363 

0.363 

13.08 

910324 

0100 

0.87 

0.79 

0.38 

14.0 

174.0 

0.103 

0.103 

0.479 

0.496 

0.496 

13.94 

910324 

0700 

0.81 

0.76 

0.28 

10.0 

-90.0 

0.113 

0.113 

0.371 

0.332 

0.332 

13.09 

910324 

1300 

0.83 

0.73 

0.37 

8.0 

176.0 

0.113 

0.113 

0.494 

0.309 

0.309 

13.73 

910324 

1900 

0.69 

0.64 

0.23 

6.0 

-90.0 

0.113 

0.113 

0.393 

0.347 

0.347 

13.05 

910323 

0100 

0.92 

0.88 

0.27 

22.0 

90.0 

0.230 

0.230 

0.309 

0.303 

0.303 

13.98 

910323 

0700 

0.90 

0.67 

0.23 

26.0 

90.0 

0.171 

0.162 

0.287 

0.302 

0.339 

13.11 

910323 

1300 

0.69 

0.66 

0.21 

10.0 

166.0 

0.103 

0.103 

0.313 

0.413 

0.413 

13.69 

910323 

1900 

0.48 

0.46 

0.13 

12.0 

162.0 

0.093 

0.103 

0.334 

0.394 

0.427 

13.09 

910326 

0100 

0.49 

0.43 

0.17 

8.0 

174.0 

0.093 

0.093 

0.383 

0.420 

0.420 

13.86 

910326 

0700 

0.43 

0.40 

0.16 

10.0 

166.0 

0.093 

0.083 

0.409 

0.393 

0.303 

13.37 

910326 

1300 

0.33 

0.30 

0.22 

6.0 

176.0 

0.093 

0.093 

0.443 

0.312 

0.312 

13.68 

910326 

1900 

0.34 

0.49 

0.22 

8.0 

-178.0 

0.093 

0.093 

0.444 

0.317 

0.317 

13.39 

910327 

0100 

0.67 

0.61 

0.26 

10.0 

174.0 

0.093 

0.103 

0.426 

0.420 

0.434 

13.63 

( Continued) 
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Data 

Tlas 

1ST 

i. 

*mo.t 

Jm 

V  v 

Xm> 

*9.1 

XP.* 

Depth 

■ 

•10327 

0700 

0.73 

0.46 

0.31 

4.0 

-178.0 

0.093 

0.093 

0.468 

0.488 

0.488 

13.47 

910327 

1300 

0.10 

0.73 

0.30 

2.0 

162.0 

0.093 

0.093 

0.407 

0.419 

0.419 

13.33 

910327 

1900 

0.73 

0.63 

0.33 

0.0 

-178.0 

0.093 

0.093 

0.498 

0.498 

0.498 

13.43 

910321 

0100 

0.63 

0.38 

0.23 

2.0 

138.0 

0.083 

0.093 

0.433 

0.407 

0.31* 

13.10 

910321 

0700 

0.6* 

0.37 

0.28 

0.0 

-17*. 0 

0.083 

0.083 

0.483 

0.499 

0.499 

13.43 

910321 

1300 

0.33 

0.31 

0.21 

-4.0 

-90.0 

0.083 

0.093 

0.400 

0.37* 

0.436 

12.95 

91032S 

1900 

0.32 

0.*6 

0.23 

•22.0 

-130.0 

0.093 

0.083 

0.307 

0.313 

0.333 

13.39 

910329 

0100 

0.  *7 

0.** 

0.18 

-s. a 

-90.0 

0.093 

0.083 

0.410 

0.403 

0.478 

12.90 

910329 

0700 

0.31 

0.  *3 

0.23 

-30.0 

-90.0 

0.093 

0.093 

0.499 

0.320 

0.320 

13.66 

910329 

1300 

0.69 

0.6* 

0.27 

-28.0 

-90.0 

0.1*2 

0.2*0 

0.431 

0.380 

0.3*6 

12.86 

910329 

1900 

0.7C 

0.67 

0.37 

-24.0 

-90.0 

0.132 

0.132 

0.3*7 

0.373 

0.373 

13.83 

910330 

0100 

o.ai 

0.73 

0.32 

-20.0 

-90.0 

0.132 

0.132 

0.42* 

0.431 

0.431 

12.8* 

910330 

0700 

0.86 

0.76 

0.40 

-22.0 

-110.0 

0.103 

0.093 

0.328 

0.3*3 

0.392 

13-9* 

910330 

1300 

1.8* 

1.76 

0.33 

32.0 

90.0 

0.171 

0.171 

0.31* 

0.298 

0.298 

11.18 

910330 

1900 

1.39 

1.30 

0.33 

22.0 

90.0 

0.132 

0.132 

0.332 

0.423 

0.423 

14.30 

910331 

0100 

1.09 

1.0* 

0.31 

10.0 

1*6.0 

0.123 

0.103 

0.300 

0.239 

0.3*0 

13.0* 

910331 

0700 

1.13 

1.09 

0.37 

12.0 

164.0 

0.123 

0.132 

0.3*3 

0.3*7 

0.433 

14.01 

910331 

1300 

1.19 

1.13 

0.29 

14.0 

90.0 

0.201 

0.201 

0.233 

0.269 

0.269 

13.01 

910331 

1900 

1.28 

1.22 

0.40 

16.0 

166.0 

0.171 

0.181 

0.326 

0.320 

0.131 

14.15 

910*01 

0100 

1.03 

1.00 

0.27 

14.0 

168.0 

0.113 

0.113 

0.271 

0.303 

0.303 

12.98 

910*01 

0700 

1.02 

0.9* 

0.39 

10.0 

172.0 

0.123 

0.123 

0.416 

0.430 

0.430 

13.85 

910*01 

1300 

0.79 

0.76 

0.23 

10.0 

160.0 

0.113 

0.093 

0.30* 

0.269 

0.3*3 

12.90 

910*01 

1900 

0.66 

0.60 

0.28 

14.0 

17*. 0 

0.093 

0.093 

0.473 

0.32* 

0.32* 

13.96 

910402 

0100 

0.6* 

0.60 

0.21 

8.0 

176.0 

0.103 

0.103 

0.3*7 

0.378 

0.378 

12.99 

910*02 

0700 

0.78 

0.73 

0.27 

8.0 

17*.  0 

0.113 

0.113 

0.377 

0.4*3 

0 . 4*3 

13.72 

910*02 

1300 

0.36 

0.33 

0.18 

8.0 

-90.0 

0.123 

0.113 

0.33* 

0.288 

0.389 

12.99 

910*02 

1900 

0.32 

0.  *9 

0.18 

6.0 

17*. 0 

0.132 

0.132 

0.362 

0.*11 

0.411 

13.83 

910*03 

0100 

0.32 

0.  *8 

0.18 

6.0 

-180.0 

0.103 

0.103 

0.37* 

0.411 

0.411 

13.1* 

910*03 

0700 

0.37 

0.33 

0.19 

4.0 

174.0 

0.113 

0.103 

0.361 

0.411 

0.435 

13.38 

910*03 

1300 

0.30 

0. 47 

0.17 

4.0 

180.0 

0.113 

0.113 

0.367 

0.386 

0.386 

11.13 

910*03 

1900 

0.33 

0.30 

0.18 

0.0 

166.0 

0.113 

0.113 

0.360 

0.379 

0.379 

13.63 

910*0* 

0100 

0.36 

0.33 

0.18 

-*.0 

160.0 

0.113 

0.083 

0.332 

0.322 

0.*01 

13.26 

910*0* 

0700 

0.61 

0.38 

0.19 

*.0 

-90.0 

0.093 

0.093 

0.336 

0.468 

0.468 

13.31 

910*0* 

1300 

0.33 

0.32 

0.18 

-* .  0 

-102.0 

0.093 

0.093 

0.3*1 

0 .  *17 

0.417 

13.18 

910*0* 

1900 

0.37 

0.34 

0.20 

-8.0 

-90.0 

0.093 

0.093 

0.36* 

0.4*2 

0.4*2 

13.  *1 

910*03 

0100 

0.33 

o.*« 

0.22 

•12.0 

•1*8.0 

0.093 

0.093 

0.4*7 

0.320 

0.320 

13.36 

910*03 

0700 

0.33 

0.  *9 

0.19 

-12.0 

-90.0 

0.093 

0.093 

0.380 

0.36* 

0.36* 

13.10 

910*03 

1300 

0.36 

0.31 

0.23 

-14.0 

-1*0.0 

0.093 

0.093 

0.46* 

0.490 

0.490 

13.20 

910*03 

1900 

0.31 

0.*7 

0.19 

-20.0 

-90.0 

0.093 

0.093 

0.406 

0.413 

0.413 

13.13 

910*00 

0100 

Q.*3 

0.40 

0.19 

-2*.0 

-1*6.0 

0.093 

0.093 

0.472 

0.496 

0.496 

13.38 

910*00 

0700 

0.** 

0.*1 

0.13 

-20.0 

-90.0 

0.093 

0.093 

0.339 

0.367 

0.367 

12.96 

910*00 

1300 

0.30 

0 .  *6 

0.21 

-12.0 

-90.0 

0.103 

0.103 

0.439 

0.310 

0.510 

13.31 

910*00 

1900 

0.30 

0.47 

0.18 

-12.0 

-100.0 

0.103 

0.103 

0.393 

0.370 

0.370 

13.09 

910*07 

0100 

0.*I 

Q.*3 

0.22 

-12.0 

-1*4.0 

0.103 

0.103 

0.30* 

0.3*9 

0.5*9 

13.32 

910*07 

0700 

o.*o 

0.37 

0.16 

-12.0 

-90.0 

0.103 

0.103 

0.42* 

0 .  *39 

0.439 

12.97 

910*07 

1300 

0.36 

0.33 

0.13 

-12.0 

-1*4.0 

0.103 

0.113 

0.468 

0.433 

0.333 

13.41 

910*07 

1900 

0.33 

0.32 

0.14 

-12.0 

•100.0 

0.113 

0.113 

0.429 

0.396 

0.396 

13.03 

(Continued) 
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(Coctlouad) 


Sam  H«l«hc  ,  IHg9«l9B  fraauancr  *a*l tat  Ion 


Data 

Tlaa 

EEL. 

*■0,1 

«-o.r 

• 

V  V  -  *>•* 

Dapth 

■ 

910*0* 

0100 

0.36 

0.33 

0.14 

-3*.0 

-10*. 0  0.113  0.113  0.431  0.46* 

0.46* 

13.37 

910*0* 

0700 

0.3* 

0.31 

0.12 

-3*.0 

-90.0  0.113  0.113  0.39*  0.373 

0.373 

12.9* 

910*0* 

1300 

0.33 

0.31 

0.1* 

•32.0 

-106.0  0.113  0.113  0.4*3  0.479 

0.479 

13.47 

910*08 

1900 

0.33 

0.32 

0.13 

•30.0 

-90.0  0.113  0.113  0.412  0.377 

0.377 

13.00 

910*09 

0100 

0.3* 

0.33 

0.16 

-s*.o 

-90.0  0.113  0.113  0.461  0.437 

0.437 

13.37 

910*09 

0700 

o.*s 

0.4* 

0.17 

-34.0 

-90.0  0.162  0.1*2  0.393  0.41* 

0.41* 

12.93 

910*09 

1300 

0.31 

0.47 

0.19 

-26.0 

-90.0  0.1*2  0.132  0.403  0.416 

0.4*2 

13.43 

910*09 

1900 

0.31 

0.4* 

0.1* 

-2*.0 

-90.0  0.113  0.113  0.377  0.3*7 

0.367 

12.93 

910*10 

0100 

0.32 

0.46 

0.23 

-26. C 

-90.0  0.1*2  0.1*2  0.4*7  0.4*9 

0.4*9 

13.4* 

910*10 

0700 

0.43 

0.42 

0.17 

-36.0 

-90.0  0.1*2  0.142  0.411  0.38* 

0.38* 

12.94 

910*10 

1300 

0.39 

0.36 

0.15 

-26.0 

-90.0  0.113  0.113  0.428  0.437 

0.437 

13.42 

O 

r4 

O 

H 

1900 

0.3* 

0.33 

0.15 

-34.0 

-90.0  0.113  0.113  0.422  0.414 

0.414 

13.17 

910*11 

0100 

0.69 

0.66 

0.22 

44.0 

90.0  0.240  0.240  0.340  0.336 

0.336 

13.61 

910*11 

0700 

l.*2 

1.7* 

0.52 

42.0 

90.0  0.162  0.171  0.299  0.21* 

0.365 

13.32 

910*11 

1300 

0.99 

0.9* 

0.29 

24.0 

162.0  0.162  0.171  0.311  0.276 

0.387 

13.59 

910*11 

1900 

0.33 

0.32 

0.19 

10.0 

164.0  0.162  0.132  0.367  0.387 

0.477 

13.38 

910*12 

0100 

0 .  *9 

0.46 

0.15 

4.0 

-172.0  0.1*2  0.171  0.331  0.33* 

0.423 

13.29 

910*12 

0700 

0.79 

0.77 

0.19 

18.0 

90.0  0.210  0.2*0  0.241  0.203 

0.233 

13.47 

910*12 

1300 

l.*0 

1.36 

0.33 

26.0 

90.0  0.162  0.162  0.2*3  0.242 

0.2*2 

13.40 

910*12 

1900 

1.03 

0.99 

0.29 

24.0 

16*. 0  0.162  0.162  0.297  0.313 

0.313 

13.66 

910*13 

0100 

1.11 

1.0* 

0.26 

14.0 

90.0  0.201  0.171  0.241  0.234 

0.307 

13.1* 

910*13 

0700 

o.a* 

o.ai 

0.23 

1*.0 

168.0  0.171  0.171  0.278  0.302 

0.302 

13.5* 

910*13 

1300 

0.80 

0.7* 

0.20 

6.0 

-9 0.0  0.162  0.162  0.236  0.2*3 

0.2*3 

13.1* 

910*13 

1900 

0.*3 

0.79 

0.26 

2.0 

-90.0  0.132  0.142  0.32*  0.2*3 

0.369 

13.81 

910*1* 

0100 

0.7* 

0.75 

0.21 

0.0 

-90.0  0.152  0.171  0.2*3  0.268 

0.293 

13.00 

910*1* 

0700 

0.63 

0.61 

0.22 

-l*.0 

-90.0  0.162  0.162  0.369  0.393 

0.395 

13.78 

910*1* 

1300 

0.33 

0.31 

0.15 

-12.0 

-90.0  0.093  0.093  0.302  0.364 

0.364 

13.16 

91041* 

1900 

1.30 

i.*l 

0.32 

22.0 

138.0  0.132  0.132  0.372  0.417 

0.417 

14.23 

910*13 

0100 

1.33 

1.30 

0.29 

14.0 

90.0  0.113  0.103  0.224  0.20* 

0.243 

13.00 

910*13 

0700 

1.28 

1.22 

0.40 

14.0 

176.0  0.113  0.113  0.323  0.365 

0.363 

13.87 

910*13 

1300 

1.00 

0.96 

0.26 

6.0 

-90.0  0.113  0.0*3  0.270  0.236 

0.2*3 

12.8* 

910*13 

1900 

o.*« 

0.77 

0.32 

0.0 

178.0  0.123  0.113  0.413  0.411 

0.4*3 

14.10 

910*1* 

0100 

0.71 

0.6* 

0.20 

-*.0 

-90.0  0.103  0.113  0.294  0.238 

0.2*3 

12.76 

910*16 

0700 

0.6* 

0.62 

0.28 

-12.0 

-90.0  0.113  0.103  0.432  0.426 

0.616 

13.84 

910*16 

1300 

0.33 

0.30 

0.16 

-10.0 

-90.0  0.123  0.123  0.327  0.334 

0.33* 

12.81 

910*16 

1900 

0.62 

0.38 

0.21 

2.0 

-90.0  0.123  0.132  0.332  0.3*6 

0.405 

14.12 

910*17 

0100 

0.37 

0.55 

0.15 

-12.0 

-90.0  0.132  0.162  0.277  0.224 

0.267 

12.75 

910*17 

0700 

0.39 

0.33 

0.20 

-12.0 

-90.0  0.162  0.123  0.356  0.333 

0.412 

13.68 

910*17 

1300 

0 ,  *7 

0.4* 

0.1* 

-1*.0 

-90.0  0.  7  0.064  0.312  0.290 

0.534 

12.73 

910417 

1900 

0 .  *7 

0.43 

0.19 

-22.0 

•90.0  0.132  0.06*  0.430  0.420 

0.7*4 

13.9* 

910*1* 

0700 

1.03 

1.00 

0.28 

18.0 

90.0  0.240  0.240  0.277  0.235 

0.233 

13.  *2 

910*1* 

1300 

1.70 

1.63 

0.42 

32.0 

90.0  0.171  0.171  0.233  0.251 

0.231 

13.37 

91041* 

1900 

1.66 

1.60 

0.  *3 

30.0 

90.0  0.152  0.132  0.2*1  0.2*0 

0.2*0 

14.11 

910*19 

0100 

1.7* 

1.73 

0.43 

22.0 

90.0  0.171  0.171  0.247  0.2*2 

0.2*2 

13.63 

910*19 

0700 

1.7* 

1.6* 

0.43 

20.0 

90.0  0.171  0.171  0.23*  0.24* 

0.24* 

13.71 

910*19 

1300 

1.33 

1.47 

0.39 

14.0 

90.0  0.103  0.103  0.267  0.290 

0.290 

13.30 

910*19 

1900 

1.39 

1.33 

0.41 

10.0 

90.0  0.103  0.103  0.2*7  0.293 

0.293 

13.72 
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1.60 

1.33 
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8.0 

-90.0 

0.1*2 
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0.238 

0.211 

0.271 

13.8* 

910*20 

0700 

1.84 

1.79 

0.42 

4.0 

-90.0 

0.1*2 

0.1*2 

0.23* 

0.223 

0.223 

13.45 

•10*20 

1300 

2.38 

2.32 

0.33 

0.0 

-90.0 

0.113 

0.113 

0.217 

0.209 

0.209 

13.77 

•10*20 

1900 

3.07 

2.99 

0.70 

12.0 

90.0 

0.123 

0.2*0 

0.23* 

0.133 

0.4*5 

13.37 

910421 

0100 

2.17 

2.08 

0.62 

16.0 

90.0 

0.123 

0.113 

0.299 

0.266 

0.313 

1*.U 

910*21 

0700 

1.66 

1.39 

0.48 

18.0 

90.0 

0.123 

0.103 

0.302 

0.230 

0.313 

13.33 

•10*21 

1300 

1.36 

1.23 

0.32 

24.0 

-178.0 

0.123 

0.123 

0.417 

0.378 

0.378 

13.92 

910*21 

1900 

0.83 

0.79 

0.32 

20.0 

-178.0 

0.123 

0.113 

0.401 

0.323 

0.394 

13.22 

910*22 

0100 

0.61 

0.33 

0.26 

18.0 

180.0 

0.123 

0.113 

0.473 

0.316 

0.523 

13.99 

910*22 

0700 

0.43 

0.41 

0.17 

6.0 

-180.0 

0.123 

0.123 

0.422 

0.423 

0.423 

13.00 

910*22 

1300 

0.46 

0.43 

0.21 

4.0 

•90.0 

0.113 

0.113 

0.48* 

0.506 

0.306 

13.73 

910*22 

1900 

0.4  9 

0.45 

0.20 

-4.0 

-90.0 

0.113 

0.113 

0.4*0 

0.422 

0.422 

12.90 

910*23 

0100 

0.37 

0.30 

0.27 

2.0 

-1*0.0 

0.113 

0.113 

0.333 

0.585 

0.385 

13.80 

910*23 

0700 

0.*7 

0.43 

0.18 

2.0 

-90.0 

0.113 

0.103 

0.422 

0.360 

0.4*5 

12.97 

910*23 

1300 

0.46 

0.41 

0.21 

6.0 

176.0 

0.123 

0.113 

0.319 

0.478 

0.558 

13.82 

910*23 

1900 

0.47 

0.43 

0.19 

4.0 

168.0 

0.123 

0.123 

0.430 

0.373 

0.373 

13.06 

910*2* 

0100 

0.31 

0.46 

0.23 

0.0 

174.0 

0.123 

0.123 

0.492 

0.452 

0.452 

13.72 

910*2* 

0700 

0.3* 

0.30 

0.21 

2.0 

-90.0 

0.123 

0.093 

0.413 

0.342 

0.4*9 

13.05 

910*23 

1300 

0.*7 

0.4* 

0.13 

6.0 

90.0 

0.2*0 

0.093 

0.3*5 

0  306 

0.476 

13.5* 

910*23 

1900 

0.43 

0-*2 

0.17 

10.0 

-180.0 

0.083 

0.093 

0.401 

0.4*0 

0.491 

13.33 

910*26 

0100 

0.  *6 

0.4* 

0.13 

10.0 

176.0 

0.093 

0.093 

0.3*3 

0.369 

0.369 

13.33 

910*26 

0700 

0.47 

0.*3 

0.18 

2.0 

-U0.0 

0.093 

0.093 

0.429 

0.498 

0.498 

13.48 

910*26 

1300 

0.30 

0.  *7 

0.18 

•2.0 

176.0 

0.1*2 

0.103 

0.382 

0.37* 

0.482 

13.32 

910*26 

1900 

0.43 

0.39 

0.18 

-2.0 

180.0 

0.093 

0.093 

0.451 

0.492 

0.492 

13.72 

910*27 

0100 

0.42 

0.39 

0.13 

-*.0 

172.0 

0.093 

0.093 

0.378 

0.515 

0.515 

13.09 

910*27 

0700 

0.43 

0.39 

0.17 

0.0 

176.0 

0.093 

0.093 

0.435 

0.4*3 

0.4*3 

13.59 

910*27 

1300 

0.3* 

0.30 

0.19 

2.0 

-90.0 

0.093 

0.093 

0.37* 

0.432 

0.  *32 

13.11 

910*27 

1900 

0.36 

0.30 

0.2* 

-2.0 

17*.  0 

0.093 

0.093 

0.470 

0.509 

0.309 

13.86 

910*28 

0100 

0.63 

0.39 

0.23 

2.0 

168.0 

0.07* 

0.07* 

0.391 

0.479 

0.479 

12.91 

910*28 

0700 

0.60 

0.33 

0.26 

0.0 

172-0 

0.083 

0.083 

0.495 

0.3*7 

0.3*7 

13.66 

910*28 

1300 

0.36 

0.32 

0.21 

0.0 

166.0 

0.083 

0.083 

0.39* 

0.371 

0.371 

12.97 

910*28 

1900 

0.33 

0.49 

0.20 

-26.0 

-90.0 

0.083 

0.083 

0.*06 

0.435 

0.433 

14.0* 

910*29 

0100 

0.32 

0.48 

0.19 

-30.0 

-90.0 

0.093 

0.083 

0.396 

0.377 

0.6*5 

13.0* 

910*29 

0700 

0.80 

0.76 

0.23 

-l*-0 

-90.0 

0.162 

0.162 

0.323 

0.3*5 

0.3*5 

13.91 

910*30 

1300 

0.6* 

0.61 

0.20 

6.0 

-90.0 

0.132 

0.07* 

0.323 

0.275 

0.417 

13.00 

910*30 

1900 

0.33 

0.49 

0.20 

*.0 

180.0 

0.07* 

0.07* 

0 .  *13 

0.472 

0.472 

14.00 

910301 

0100 

0.46 

0.42 

0.18 

2.0 

168.0 

0.07* 

0.07* 

0.427 

0.591 

0.591 

13.00 

910301 

0700 

0.41 

0.37 

0.16 

4.0 

172.0 

0.07* 

0.07* 

0.435 

0.603 

0.603 

13.7* 

910301 

1300 

0.40 

0.37 

0.16 

8.0 

168.0 

0.07* 

0.07* 

0.444 

0.379 

0.379 

13.08 

910301 

1900 

0.38 

0.33 

0.13 

-4.0 

180-0 

0.083 

0.083 

0.*1* 

0.5*6 

0.5*6 

13.96 

910302 

0100 

0.33 

0.29 

0.1* 

-2.0 

176.0 

0.083 

0.083 

0.492 

0.335 

0.333 

13.13 

910302 

0700 

0.33 

0.31 

0.12 

2.0 

172.0 

0.083 

0.07* 

0.39* 

0.386 

0.506 

13.65 

910302 

1300 

0.3* 

0.31 

0.1* 

-2.0 

172.0 

0.083 

0.07* 

0.436 

0.*«6 

0.588 

13.15 

910302 

1900 

0.37 

0.3* 

0.13 

-4.0 

172.0 

0.07* 

0.07* 

0.390 

0.493 

0.493 

13.80 

910303 

0100 

0.38 

0.33 

0.1* 

-2.0 

176.0 

0.132 

0.07* 

0.401 

0.298 

0.301 

13.2* 

910303 

0700 

0.31 

0.29 

0.11 

0.0 

172.0 

0.083 

0.07* 

0.39* 

0.401 

0.62* 

13.  *7 
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■ 

•10203 

1300 

0.«2 

0.39 

0.14 

18.0 

90.0 

0.240 

0.240 

0.346 

0.328 

0.328 

13.30 

910303 

1900 

0.37 

0.55 

0.13 

46.0 

90.0 

0.191 

0.201 

0.270 

0.238 

0.263 

13.66 

•1030* 

0100 

0.41 

0.40 

0.11 

36.0 

166.0 

0.181 

0.181 

0.288 

0.236 

0.236 

13.39 

•1030* 

0700 

0.30 

0.28 

0.10 

4.0 

-90.0 

0.240 

0.083 

0.332 

0.267 

0.486 

13.36 

91030* 

1300 

0.30 

0.29 

0.09 

-6.0 

-90.0 

0.191 

0.191 

0.323 

0.367 

0.367 

13.37 

910304 

1900 

0.30 

0.29 

0.09 

-2.0 

-90.0 

0.201 

0.093 

0.304 

0.261 

0.523 

13.43 

910303 

0100 

0.40 

0.39 

0.11 

14.0 

-90.0 

0.171 

0.171 

0.286 

0.284 

0.284 

13.47 

910303 

0700 

0.40 

0.39 

0.10 

18.0 

-90.0 

0.191 

0.191 

0.270 

0.249 

0.249 

13.20 

910303 

1300 

0.37 

0.35 

0.11 

-14.0 

-90.0 

0.230 

0.083 

0.310 

0.274 

0.378 

13.48 

910303 

1900 

0.41 

0.39 

0.12 

-6.0 

-90.0 

0.191 

0.191 

0.299 

0.307 

0.307 

13.39 

910306 

0100 

0.49 

0.46 

0.16 

-18.0 

-90.0 

0.201 

0.201 

0.332 

0.383 

0.383 

13.69 

910306 

0700 

0.68 

0.64 

0.21 

-18.0 

-90.0 

0.181 

0.181 

0.319 

0.333 

0.333 

13.17 

910306 

1900 

0.39 

0.38 

0.12 

-14.0 

-90.0 

0.171 

0.132 

0.317 

0.248 

0.372 

13.15 

910307 

0100 

0.40 

0.37 

0.14 

-22.0 

-138.0 

0.162 

0.152 

0.386 

0.343 

0.381 

13.62 

910507 

0700 

0.33 

0.33 

0.12 

-20.0 

-90.0 

0.152 

0.142 

0.361 

0.320 

0.348 

13.11 

910507 

1300 

0.68 

0.66 

0.18 

62.0 

90.0 

0.240 

0.240 

0.278 

0.262 

0.262 

13.77 

910307 

1900 

0.50 

0.48 

0.14 

28.0 

90.0 

0.210 

0.220 

0.283 

0.239 

0.330 

13.17 

910506 

0100 

0.45 

0.43 

0.14 

0.0 

-90.0 

0.210 

0.083 

0.336 

0.211 

0.337 

13.68 

910509 

0700 

0.41 

0.39 

0.12 

0.0 

130.0 

0.142 

0.093 

0.303 

0.246 

0.351 

13.00 

910308 

1300 

0.40 

0.38 

0.14 

0.0 

180.0 

0.093 

0.093 

0.358 

0.435 

0.435 

13.71 

910308 

1900 

0.40 

0.39 

0.13 

-8.0 

-90.0 

0.162 

0.093 

0.323 

0.282 

0.412 

13.08 

910509 

0100 

0.42 

0.40 

0.1* 

-10.0 

-144.0 

0.181 

0.171 

0.363 

0.339 

0.410 

13.66 

910309 

0700 

0.42 

0.40 

0.12 

-14.0 

-90.0 

0.191 

0.093 

0.298 

0.230 

0.460 

12.98 

910509 

1300 

0.38 

0.36 

0.12 

-14.0 

174.0 

0.191 

0.093 

0.336 

0.296 

0.485 

13.73 

910509 

1900 

0.40 

0.38 

0.10 

>12.0 

-90.0 

0.103 

0.103 

0.266 

0.288 

0.288 

13.18 

910310 

0100 

0.47 

0.4* 

0.14 

-10.0 

-90.0 

0.201 

0.201 

0.312 

0.290 

0.290 

13.66 

910310 

0700 

0.47 

0.43 

0.13 

-18.0 

-90.0 

0.201 

0.191 

0.293 

0.23* 

0.326 

13.08 

910510 

1300 

0.45 

0.41 

0.16 

-14.0 

-1*4.0 

0.132 

0.132 

0.397 

0.491 

0.491 

13.67 

910510 

1900 

0.39 

0.37 

0.13 

-18.0 

-142.0 

0.132 

0.132 

0.408 

0.480 

0.480 

13.33 

910311 

0100 

0.93 

0.91 

0.21 

10.0 

90.0 

0.220 

0.201 

0.230 

0.193 

0.289 

13.33 

910511 

0700 

0.93 

0.93 

0.20 

10.0 

90.0 

0.191 

0.1*1 

0.217 

0.206 

0.206 

13.28 

910511 

1300 

0.93 

0.89 

0.23 

8.0 

146.0 

0.123 

0.123 

0.284 

0.303 

0.305 

13.37 

910511 

1900 

0.69 

0.66 

0.20 

2.0 

172.0 

0.123 

0.123 

0.302 

0.329 

0.329 

13.51 

910312 

0100 

0.7* 

0.71 

0.22 

2.0 

166.0 

0.103 

0.113 

0.303 

0.343 

0.345 

13.20 

910312 

0700 

0.73 

0.68 

0.28 

-4.0 

•148.0 

0.113 

0.113 

0.412 

0.413 

0.413 

13.31 

910512 

1300 

0.78 

0.73 

0.27 

-6.0 

•146.0 

0.103 

0.103 

0.377 

0.409 

0.409 

13.23 

910312 

1900 

0.62 

0.36 

0.27 

•6.0 

-148.0 

0.103 

0.103 

0.474 

0.463 

0.463 

13.67 

910313 

0100 

0.39 

0.36 

0.19 

-6.0 

162.0 

0.103 

0.103 

0.337 

0.323 

0.323 

12.90 

910313 

0700 

0.32 

0.47 

0.22 

-4.0 

178.0 

0.093 

0.103 

0.471 

0.440 

0.479 

13.46 

910313 

1300 

0.27 

0.23 

0.09 

•42.0 

-90.0 

0.171 

0.103 

0.363 

0.300 

0.557 

12.76 

910313 

1900 

0.29 

0.27 

0.09 

-52.0 

-90.0 

0.181 

0.191 

0.329 

0.336 

0.346 

14.27 

910316 

0100 

0.37 

0.36 

0.10 

-38.0 

•90.0 

0.240 

0.132 

0.277 

0.219 

0.333 

12.92 

910316 

0700 

0.34 

0.51 

0.18 

20.0 

90.0 

0.220 

0.210 

0.332 

0.309 

0.336 

13.88 

910316 

1300 

0.48 

0.47 

0.12 

10.0 

90.0 

0.230 

0.230 

0.262 

0.231 

0.251 

12.83 

910316 

1900 

0.39 

0.37 

0.17 

-6.0 

176.0 

0.171 

0.162 

0.306 

0.332 

0.342 

14.13 

910317 

0100 

0.46 

0.44 

0.11 

-6.0  -90.0  0.162 
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■ 

*10317 

0700 

0.** 

0.42 

0.13  -4.0  -90.0 

0.201 

0.210 

0.301 

0.267 

0.304 

13.70 

*10317 

1300 

0.36 

0.33 

0.10  -12.0  -*0.0 

0.171 

0.103 

0.213 

0.211 

0.435 

12.97 

*10317 

1900 

0.30 

0.36 

0.12  -6.0  -90.0 

0.191 

0.111 

0.335 

0.219 

0.326 

13.80 

•10310  0100 

0.32 

0.30 

0.10  -12.0  -90.0 

0.113 

0.113 

0.343 

0.3*7 

0.3*7 

13.19 

•10310  0700 

0.31 

0.29 

0.11  0.0  -90.0 

0.113 

0.123 

0.364 

0.394 

0.429 

13.40 

•10310 

1300 

l.*0 

1.37 

0.21  11.0  90.0 

0.201 

0.201 

0.203 

0.236 

0.236 

13.42 

•10310 

1900 

2.21 

2.16 

0.46  22.0  90.0 

0.142 

0.142 

0.213 

0.186 

0.186 

13.87 

910319 

0100 

2.22 

2.17 

0.41  20.0  90.0 

0.113 

0.113 

0.223 

0.243 

0.243 

13.93 

910319 

0700 

2.*« 

2.** 

0.42  12.0  90.0 

0.132 

0.171 

0.171 

0.150 

0.110 

13.63 

910319 

1300 

2.31 

2.26 

0.41  14.0  90.0 

0.142 

0.142 

0.210 

0.197 

0.197 

13.92 

910319 

1900 

1.17 

1.13 

0.36  12.0  90.0 

0.162 

0.171 

0.199 

0.191 

0.207 

13.33 

910320 

0100 

1.63 

1.60 

0.39  6.0  176.0 

0.093 

0.093 

0.242 

0.199 

0.199 

14.04 

910320 

0700 

1.3* 

1.31 

0.32  2.0  -90.0 

0.093 

0.171 

0.210 

0.179 

0.243 

13.30 

910321 

1300 

1.00 

0.93 

0.29  -4.0  -90.0 

0.123 

0.132 

0.310 

0.279 

0.325 

13.84 

910321 

1900 

0.92 

O.M 

0.26  -4.0  -90.0 

0.123 

0.162 

0.294 

0.230 

0.352 

12.97 

910322 

0100 

0.16 

0.10 

0.31  -6.0  -90.0 

0.132 

0.132 

0.383 

0.392 

0.392 

13.76 

910322 

0700 

0.73 

0.72 

0.22  -1.0  -90.0 

0.123 

0.142 

0.307 

0.256 

0.311 

12.8* 

910322 

1300 

0.67 

0.62 

0.23  -1.0  -1*4.0 

0.132 

0.113 

0.400 

0.392 

0.412 

13.79 

910322 

1900 

0.37 

0.34 

0.20  -1.0  -90.0 

0.113 

0.123 

0.371 

0.323 

0.360 

12.96 

910323 

0100 

0.37 

0.32 

0.24  -6.0  -144.0 

0.142 

0.123 

0.469 

0.463 

0.496 

13.39 

910323 

0700 

0.  *7 

0.*3 

0.17  -20.0  -140.0 

0.123 

0.113 

0.402 

0.326 

0.389 

12.91 

910323 

1300 

o.*a 

0.43 

0.22  -26.0  -144.0 

0.123 

0.123 

0.511 

0.567 

0-567 

13.76 

910323 

1900 

o.*o 

0.36 

0.17  -24.0  -131.0 

0.123 

0.113 

0.479 

0.442 

0.516 

13.24 

91032* 

0100 

o.*o 

0.36 

0.17  -21.0  -106.0 

0.123 

0.123 

0.453 

0.4*5 

0.443 

13.62 

91032* 

0700 

0.  *2 
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Figure  Bl.  Bulk  dace  for  November  1990 
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Figure  B2.  Bulk  data  for  Deceober  1990 


BA 


Januory,  1991 

1  2  3  4  5  6  7  8  9  10  IT  12  13  14  IS  t6  17  IB  19  20  21  22  2  3  24  2  3  26  27  28  2  9  3  0  31 


Figura  B3.  Bulk  data  for  January  1991 
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Figure  B4.  Bulk  data  for  February  1991 
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Figure  B5.  Bulk  data  for  March  1991 
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Figure  B6.  Bulk  data  for  April  1991 
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Figure  B7.  Bulk  data.  for  Hay  1991 
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Figure  B8.  Bulk  dace  for  June  1991 
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Appendix  C:  A  Derivation  of  the  Maximum  Likelihood  Estimator 


1.  Under  ideal  circumstances,  a  single  train  of  ocean  surface  waves  can 
be  represented  as  a  sinusoid  in  the  form* 

ij(x,7,t)  ■  a  cos(kcos&  x  *  k sintfy  -  ot  ♦  «)  (Cl) 

where  x  and  y  are  ixght-handed  coordinates  in  a  horizontal  plane  (at 
elevation  z  -  0,  where  z  is  the  vertical  coordinate  directed  positive 
upward  from  the  mean  sea  surface),  t  is  time,  a  is  the  wave  amplitude,  k 
is  radian  wavenumber,  8  is  the  direction  toward  which  the  wave  propagates 
and  is  measured  counterclockwise  (looking  downward)  from  the  positive  x-axis, 
a  is  radian  wave  frequency,  and  e  is  an  initial  phase.  When  linear  wave 
theory  is  used,  wavenumber  and  frequency  are,  for  a  given  depth  of  water, 
uniquely  related  by  the  linear  wave  dispersion  relation 

a2  ■  gktanhkd  (C2) 

where  g  is  gravitational  acceleration,  and  d  is  water  depth.  A  natural 
sea  can  be  considered  an  infinite  sum  of  such  waves  so  that  sea  surface 
displacement  is  represented  generally  as 


if (x,y,t)  -  £  aacos(khcos8nx  *  *^sin*ny  -  ont  ♦  «n) 

n«l 


(C3) 


For  each  index  n  ,  there  is  a  unique  set  of  a„  ,  trn  ,  8a  and  <„  that,  for 
a  given  depth,  completely  determine  the  nature  of  the  n**  wave  train  (note 
that  ICq  is  dictated  by  an  and  d  through  the  dispersion  relation) . 

2.  because  it  can  be  rather  awkward  keeping  track  of  this  infinite 
number  of  four-variable  sets,  it  is  useful  to  have  a  form  of  book-keeping  that 
is  compact  enough  that  even  a  very  complicated  sea  can  be  quantified  and 
depicted  in  a  simple  way.  A  common  form  of  such  book-keeping  is  the 
frequency- direction  spectrum.  In  this  method,  a  plane  is  gridded  in  even 
steps  of  the  orthogonal  coordinates  a  and  i  ,  with  each  cell  in  the  grid 
being  of  infinitesimal  dimension  to  by  LB  .  There  is  then  a  cell  that 
corresponds  to  the  an  and  8n  of  the  n*1*  wave  train  so  that  the  other  two 
properties  a,,  and  <n  of  this  wave  train  can  be  associated  with  that  cell. 
Above  the  gridded  a-8  plane  can  be  constructed  surfaces  that  represent  wave 


A  notation  for  this  appendix  is  included  in  Appendix  D 
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Crain  properties.  The  mosc  common  of  Chese  is  variance  of  sea  surface 
displacement  (which  by  linear  wave  theory  is  proportional  to  mean  energy  of  a 
wave  train)  normalized  by  grid  cell  area.  This  variance  density  (sometimes 
called  energy  density)  together  with  its  coordinates  in  o-S  space  is  the 
frequency -direction  spectrum  of  sea  surface  variance  E(<7,0)  ,  which  is  most 
commonly  sought.  One  could  also  build  spectra  of  wave  amplitude  or  phase,  but 
wave  energy  is  usually  of  greatest  interest. 

3.  The  variance  of  a  single  sinusoidal  wave  train  of  amplitude  a„  is 
half  the  amplitude  squared.  By  definition,  this  is  the  volume  under  the 
spectrum  over  one  grid  cell.  This  relation  can  be  written 


K-j 


(C4) 


The  primary  effort  in  most  directional  spectral  estimation  is  the  determina¬ 
tion  of  E (.0,6)  from  time  series  of  sea  surface  displacement  at  a  few  points 
in  space,  i.e.,  from  a  sparse  sample  of  the  left  side  of  Equation  C3  to  an 
estimate  of  the  integrand  on  the  right  side  of  Equation  C4. 

4.  Preliminary  theoretical  considerations  involve  computation  of  a 
Fourier  frequency  transform  of  the  covariance  between  the  model  sea  surface  at 
location  x,y  and  time  t  and  the  sea  surface  at  arbitrarily  displaced 
coordinates  x+X,y+Y  and  displaced  time  t+T  .  This  process  yields  a  useful 
relationship  between  the  measurable  frequency  cross-spectra  between  two  points 
in  space  and  an  integral  with  respect  to  direction  of  the  frequency-direction 
spectrum.  A  good  discussion  of  this  process  is  given  by  Munk,  et  al. 

(1963)*,  and  the  following  is  based  on  insights  gained  therefrom. 

5.  The  covariance  R(X,Y,T)  between  spatially  and  temporally  separated 
sea  surface  displacements  is  found  from  the  average  over  large  reaches  of 
space  and  long  durations  the  product  of  the  model  for  »j(x,y,t)  with  the 
model  for  rj(x+X,y+Y, t+T)  .  This  can  be  expressed  as 

11m 

R(X,Y,T)  -  ^  _L  J\  J_\  J*  rj(x+X  ,y+Y , e+T)  dtdxdy  («) 

C-** 


When  the  model  Equation  C3  is  substituted,  this  expression  becomes  a  triple 


References  for  this  appendix  are  listed  in  the  main  body  of  the  report 
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Integral  of  the  product  of  two  Infinite  subs.  Its  evaluation  is  time  consum¬ 
ing  but  straightforward,  and  results  in  the  rather  simple  expression 

R(X,Y,T)  -  £  14  cos<J^cos*BX  *  ^sintfnf  -  aaT )  <C6) 

a-l  2 


The  variance  within  the  sum  can  be  replaced  with  its  spectral  definition  given 
by  Equation  C4,  which  leads  to 


R(.X,Y,T ) 


I 

n*l 


cos(Jc^cos0nX  ♦  l^sin0nY  -  oaT)  (C7) 


The  summation  covers  all  conceivable  bins  in  the  infinitesimally  gridded  <7-0 
domain  and  effectively  links  the  edges  of  each  integration  region  to  form  a 
continuous  double  integral  over  the  bounds  (0  <  8  £  2*  ,  0  5  a  s  »)  of 
frequency-direction  space.  Hence,  the  covariance  function  becomes 

R(X,Y,T )  "  £*  £  £(<7*  *)  cos  (kcosSX  *  kslnB  Y  -  oT)  da  dB  (C8) 


The  trigonometric  identity  for  the  cosine  of  the  sum  of  two  angles  can  be  used 
to  isolate  terms  in  oT  in  Equation  C8.  If  that  result  Is  regrouped  and  the 
order  of  integration  reversed,  the  result  is 


-nr 


£(<7,0)  cos(kcos0  X  +  ksinff  Y)  dB 


cos  oT  da 


♦jr[r  £(<7,0 )  sin(kcos0  X  +  ksinB  Y)  dB  1  sin oT  do 


(C9) 


The  bracketed  term  in  the  first  integral  on  the  right  side  of  Equation  C9  is 
the  usual  definition  of  the  coincident  spectrum,  denoted  C(X,Y ,a)  ,  and  the 
bracketed  term  in  the  second  integral  on  the  right  side  of  Equation  C9  is  the 
quadrature  spectrum  Q(X,Y,o)  .  The  covariance  function  then  becomes 

R(X,Y,T)  -  £  [C(X,Y,a)  cos  oT  ♦  Q(X,Y,a)  sinaT]  da  (CIO) 

A  Fourier  transform  of  the  covariance  function  acts  to  isolate  the  coincident 
anci  quadrature  spectra  in  the  form 

£j?(X,Y,T)  e*UIdr  •  C(X,Y,<r)  -  lQ(X,Y,o)  (CU) 

If  the  integral  definitions  of  these  spectra  are  substituted,  one  obtains 
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C{X,Y,a)  -  LQ{X,Y,a) 


•  [*  E(a,B)  cos(kcos6  X  ♦  kstnB  Y)d6  -  i[**  E(a,9)  sin(fccos#  X  *  ksinB  Y)  dB 
Jo  Jo  (C12) 

■  E(a,8)  [  cos(kcos$  X  *  ksind  Y)  -  l  sin(fccos0X  ♦  ksinB  7)]  dB 

-  E(o,e)e-iace°“x +  dS 


6.  The  cross-spectra  can  usually  be  estimated  from  measurements  because 
samples  can  be  obtained  often  enough  that  frequency  content  is  well -resolved. 
If  samples  could  be  obtained  in  space  with  equivalent  coverage,  Equation  C12 
could  be  Fourier  transformed  spatially  to  isolate  the  frequency-direction 
spectrum.  The  problem  is  that  it  is  very  difficult  to  obtain  spatial  sampling 
with  sufficient  density  for  Fourier  transformation  and  only  a  few  (order  10) 
spatial  points  are  measured.  An  alternative  approach  is  to  create  a  model  of 
E (a, 8)  and  employ  the  sparse  available  measurements  to  determine  the  model 
parameters.  There  are  several  ways  to  formulate  such  a  model.  Maximum 
Likelihood  Estimation  (MLE)  is  one  of  them. 

7.  In  MLE,  the  combined  amplitude  and  phase  of  a  single  wave  train  are 
modeled.  This  wave  train  has  a  unique  frequency  and  direction,  and  so  has  a 
unique  location  in  the  frequency-direction  plane  used  to  define  spectral 
elements.  Spectral  density  is  computed  from  the  amplitude  of  this  wave  train. 
MLE  Is  somewhat  unique  in  that  continuous  distributions  of  energy  density  are 
not  found.  The  estimator  "looks"  along  a  single  discrete  direction  at  a  time 
to  find,  for  each  of  a  set  of  discrete  frequencies,  the  most  likely  (maximum 
likelihood)  amplitude  of  a  wave  moving  in  that  direction.  By  looking  along  a 
number  of  adjacent  discrete  directions,  an  estimated  directional  spectrum  can 
be  constructed.  Some  features  of  this  method,  including  its  direction 
resolving  ability,  its  ability  to  distinguish  multiple  wave  trains  and  its 
behavior  in  relation  to  the  number  and  spacing  of  the  observing  gages,  must  be 
deduced  empirically  by  analyzing  synthetic  data  imposed  on  a  particular  array 
design. 

8.  The  purpose  here  is  to  outline  the  reasoning  and  mathematical 
underpinnings  of  MLE  as  applied  to  wave  directional  estimation.  Vhat  follows 
is  an  expansion  of  what  is  covered  in  Just  a  few  lines  in  each  of  two  papers 
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used  for  reference.  One  is  by  Davis  and  Regier  (1977)  and  Che  ocher  is  by 
Oltman-Shay  and  Guza  (1984). 

9.  The  frequency- direccion  spectrum  is  found  by  crying  Co  isolace  che 
amplitudes  of  single  plane  waves,  one  wave  crain  ac  a  cine.  A  single  plane 
wave  can  be  defined  by  Equacion  Cl,  buC  for  Che  presenc  mechod,  ic  is  desired 
Co  rewrice  this  expression  using  che  Euler  nocacion  definition  of  cosine 
[cos  x  -  (e1*  +  e~iz)/2]  7  that  Equation  Cl  becomes 


ti(x,r,t) 


_l|  ei(kcos0s  ♦  k*in«y  -  ft  *  c)  +  g-i(kcoa6x  *  kiintf y  -  <rt  •  c) 
2L 


— a  eif  ♦  k*in« y  -  ft)  +  g-i«  g-Mkcoafz  *  kslaSy  -  ft)  (C13) 

2  1 

1  d  glOteo »flx  *  k«in8y  -  ft)  +  1^»  g-i(keoi6x  -  kain0y  •  <rt) 

2  2 


where  D  -  ae1*  is  known  as  the  complex  amplitude  and  the  asterisk  means 
complex  conjugate.  Because  this  is  only  a  change  in  notation  and  not  in 
substance,  wave  train  properties  given  above  will  remain  the  same,  but  the 
notation  will  change.  Where  a  squared  amplitude  appeared  (as  in  covariance 
estimates)  previously  as  a2  ,  ic  will  now  be  DD*  .  The  infinite  suite  of 
wave  trains  that  defines  a  sea  surface  has  the  notation 


V (x,7.t) 


i<*nco»M  ♦  kn*infiny  -  arnt> 


f  1E*e'i<knCO**n*  *  kn,in*i>Jr  • 
1  ° 


(C 14) 


Then  the  variance  of  the  n“  component  wave  is  from  Equation  C4 

-Sq2  ■  -DnDn 
2  2  n  D 


(C15) 


and  Che  covariance  function  R(X,Y,T)  is 


R(X,Y,T)  -  £  1d£  cosU^cos^X  +  k^sin^Y  -  oj) 
0-1  * 


(C16) 


All  subsequent  expressions,  including  Chose  for  Che  conCinuous  frequency- 
direction  spectrum  and  Che  cross-spectra,  will  be  identical  to  those  given 
above . 

10.  Data  are  expected  to  be  some  measure  of  sea  surface  displacement  rj 
which  may  be  direct  measurements  or  derived  from  surface  corrected  time  series 
of  pressure  recorded  at  some  level  in  the  water  column.  Here  it  is  assumed 
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that  direct  measurements  have  been  made.  There  will  be  in  general  a  set  of  N 

gages,  indexed  n  -  1,2 . N  ,  from  each  of  which  is  collected  an  R-point 

time  series  sampled  at  time  interval  At  for  a  duration  of  RAt  .  Gages  are 
deployed  in  a  two-dimensional  array  in  a  pattern  that  follows,  in  two  dimen¬ 
sions,  the  guidance  given  quite  nicely  by  Davis  and  Regier  (1977)  for  one 
dimensional  arrays.  The  frequency  content  of  the  data  is  needed,  so  a 
conventional  discrete  Fourier  transform  is  computed.  For  the  n1*  gage, 
located  at  horizontal  position  x^yn  ,  the  Fourier  transform  at  the  q1* 
frequency  is,  in  accordance  with  procedures  given  by  Bendat  and  Piersol  (1971) 
or  Jenkins  and  Watts  (1968), 

VWn>  "  £  e'i2"'ir/R  q- 1,2,*  (C17) 

r*l  ^ 


where  qr  is  the  data  value  at  the  r*1*  time  step.  The  model  of  wave  ampli¬ 
tude  D  is  based  on  the  coefficients  of  the  Fourier  series  representation  of 
sea  surface  displacement.  Fourier  series  coefficients  can  be  found  directly 
from  Fourier  transforms.  For  the  q“*  frequency  aq  of  the  transformed  data 
at  location  xn,yn  ,  the  Fourier  coefficient  can  be  written 

VWa)  -  "  |VWn)  (C18 


which  is  valid  at  all  frequencies  except  zero  and  the  Nyquist,  which  two  cases 
are  ignored  with  very  little  penalty  in  this  derivation. 

11.  The  heart  of  the  maximum  likelihood  method  is  finding  an  expression 
to  represent  the  complex  amplitude  D  -  aeu  of  a  single  wave  train  at 
frequency  aq  (as  dictated  in  practice  by  the  Fourier  transform  above)  and 
traveling  in  the  single  look  direction  a  .  Because  Fourier  coefficients 
intrinsically  contain  amplitude  and  phase  information  on  a  frequency-by¬ 
frequency  basis,  it  is  reasonable  to  formulate  a  model  for  ae1'  as  a  linear 
combination  of  the  N  Fourier  series  coefficients  at  frequency  aq  deter¬ 
mined  from  the  N  gages  in  the  array.  A  linear  combination  means  a  weighted 
sum,  so  the  model  cakes  the  form 
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where  Che  wn(a)  are  a  sec  of  N  coop  lex  veighCs  unique  Co  che  look  direc- 
Cion  a  ,  which  will  be  determined  by  Che  optimization  scheme  described  below. 
The  variance  of  che  modeled  wave  Crain  is  given  by  Equacion  C15,  and,  if  che 
model  Equacion  C19  is  subsciCuCed,  Che  variance  becomes 
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where  che  summing  index  is  changed  on  che  second  sum  Co  discinguish  che 
various  Cerms  after  Che  multiplicacion  is  carried  out.  Taking  Che  complex 
conjugate  of  Che  second  sum,  carrying  out  Che  Cerm-by-term  multiplication  and 
recombining  che  results  leads  to 
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The  contribution  of  this  variance  to  the  frequency-direction  spectrum  is 
unique  in  that  it  represents  Chat  piece  of  the  total  variance  in  a  region  of 
area  A 8  by  A a  centered  at  direction  o  and  frequency  aq  in  the 
frequency- direction  plane.  When  normalized  by  the  incremental  area  AoAtf  of 
the  frequency -direction  plane,  an  estimate  of  che  frequency-direction  spectrum 
£(aq,a)  due  to  this  wave  train  can  be  written 
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where  the  hat  (A)  means  an  estimate,  and  terms  in  che  last  equality  have 
simply  been  regrouped. 
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12.  The  br&ckeced  term  in  Equation  C22  is  a  conventional  discrete 
frequency  cross -spectrum,  denoted  £^(<7,,)  and  identically  equal  to 

-  iQ„(ffq)  ,  as  described  by  Bendat  and  Piersol  (1971).  Within  the 
brackets,  the  first  Fourier  transform  term  ^(x^y,,)  is  from  a  time  series 
collected  from  gage  n  of  the  array.  The  second  Fourier  transform  term 

the  conjugate  of  the  transform  of  a  time  series  collected  from 
gage  m  ,  which  is  generally  different  from  gage  n  .  The  frequency  increment 
La  is  known  from  conventional  discrete  Fourier  transformation  to  be  La  - 
2ir/RAt  .  With  these  substitutions,  the  bracketed  term  in  Equation  C22  becomes 

R2La  H  '  (C23) 

-  £„(*,)  -  !<?„(<*(,) 

and  the  estimate  £(oq,a)  becomes 

-  ,  »  *  .  . 

£(cr  a)  >_  J  £  wn(a)  wB*(a)  £^(<7  )  -  iQ^COJ  (C24) 

13.  It  is  important  to  note  at  this  point  that  a  variance  density 
estimate  has  been  formulated  from  a  basic  Fourier  transform  definition  and  a 
more  reliable  estimate  must  be  formed  to  avoid  arithmetic  problems  later. 

Cross -spectral  estimates  formed  from  raw  Fourier  transforms  have  broad 
confidence  intervals,  i.e.,  alone  they  have  rather  high  uncertainty.  Confi¬ 
dence  is  improved  by  obtaining  long  time  series,  segmenting  records  before 
transformation,  computing  cross-spectral  estimates  for  corresponding  pairs  of 
segments,  and  averaging  segment  estimates  together.  Further  reliability  is 
obtained  when  estimates  from  adjacent  frequency  bands  are  also  averaged 
together.  The  trade-off  is  that  the  final  resolution  frequency  bandwidth  is 
broader  than  for  the  raw  estimates.  However,  sampling  strategies  can  usually 
keep  the  new  bandwidth  tolerably  narrow.  These  smoothing  methods  are  standard 
procedures  in  spectral  estimation.  The  point  here  is  that  it  is  essential  for 
smooched  estimates  to  be  used  in  maximum  likelihood  estimation. 

14.  In  this  derivation,  smoothing  can  take  place  as  soon  as  products  of 
Fourier  transforms  are  formed.  Hence,  one  can  smooth  any  of  Equations  C21  to 
C24,  under  the  assumption  that  the  weights  w(a)  are  not  affected  by  this 
process.  If  only  segment  averaging  is  done,  the  meaning  of  the  notation  in 
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the  rest  of  this  derivation  is  unchanged.  If  band  averaging  is  done,  subsums 

of  terms  will  be  computed  over  small  ranges  of  the  subscript  q  and  the 

normalizing  bandwidth  La  will  be  larger  (equal  to  the  number  of  bands 
averaged  times  the  raw  bandwidth  2*/Rdt) .  If,  for  the  remainder  of  this 
derivation,  subscript  q  is  considered  to  be  the  index  of  this  new  set  of 
frequency  bands,  then  the  notation  need  not  be  changed. 

15.  A  logical  next  step  is  to  define  LB  ,  but  this  is  difficult  to  do 

at  this  stage.  Hypothetically,  the  method  seeks  a  single  wave  train  having  a 
single  propagation  direction  a  .  Under  ideal  conditions,  the  arc  of  direc¬ 
tion  in  the  frequency-direction  domain  represented  by  this  wave  train  is  of 
zero  width,  i.e.,  the  spectrum  is  a  spike  or  delta  function  at  direction  a  . 
However,  directional  resolving  ability  is  limited  by  the  number  of  observa¬ 
tions  available,  so  there  is  a  finite  arc  width  associated  with  the  fitted 

wave  train,  just  as  there  is  a  finite  bandwidth  in  the  frequency  domain. 

Because  of  the  intricate  arithmetic  and  dependence  on  the  number  and  separa¬ 
tion  of  gages,  it  is  usually  not  possible  to  assign  a  resolution  arc  width 
until  the  estimator  is  completely  derived.  At  that  time,  a  resolution  arc 
width  can  be  found  empirically  by  examining  the  minimum  angular  separation  of 
two  wave  trains  that  can,  by  some  criterion,  be  distinguished  with  the 
estimator.  This  procedure  is  usually  done  with  artificial  wave  signals, 
solving  the  estimator  equations  for  a  large  number  of  individual  look  direc¬ 
tions  (a's)  and  then  deciding  how  small  an  arc  is  required  for  estimates  from 
adjacent  look  directions  to  be  reasonably  statistically  independent.  Results 
will  be  a  function  of  frequency  (long  waves  are  more  difficult  to  resolve)  and 
may  be  a  function  of  direction,  depending  on  the  geometry  of  the  gage  array. 

16.  Given  that,  for  now,  LB  is  not  defined,  the  working  entity  is  the 
frequency-direction  spectrum  times  an  unspecified  arc  width.  Hence,  multiply¬ 
ing  both  sides  of  Equation  C24  by  LB  results  in  what  is  essentially  the  q^ 
element  of  a  slice  of  a  direction-bounded  frequency  spectrum  of  arc  width  LB 
centered  along  look  direction  a  ,  or 

n  h 
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The  model  Is  a  weighted  sum  of  the  frequency  cross-spectra  of  sea  surface 
displacement  at  frequency  aq  between  all  possible  pairs  of  gages  n  and 
m  ,  where  the  weights  are  unique  to  the  single  direction  a  .  Weights  are 
determined  based  on  reasoning  derived  from  consideration  of  ideal  spectra. 

17.  Complete,  ideal  cross-spectra  are  related  to  the  true  frequency- 
direction  spectrum  E(<7,0)  (no  hat)  by  Equation  C12.  If  this  expression  is 
adapted  to  an  array  of  gages  positioned  arbitrarily  (for  now)  in  x,y 
coordinates,  separations  between  gages  m  and  n  are 


Xm*»-Xn 
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Use  of  Equations  C26  in  Equation  C12  gives  an  expression  for  cross -spectra  of 
sea  surface  displacement  in  terms  of  gage  separation 


.  IT  E(arB)  d6 
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18.  The  next  step  is  to  isolate  that  part  of  the  spectrum  near  look 
direction  a  .  This  is  done  by  asserting  that  the  total  directional  spectrum 
at  frequency  aq  is  made  up  of  two  parts:  a  spike  or  delta  function  at  the 
look  direction  a  that  represents  a  single,  signal-containing  wave  train,  and 
a  "noise"  spectrum  that  represents  some  true  noise,  but  also  the  rest  of  the 
true  spectrum  that  is  to  have  minimum  influence  on  the  estimate  at  look 
direction  a  .  That  is,  a  kind  of  filter  is  enforced  that  only  allows 
information  near  look  direction  o  .  The  two  parts  of  E(aq,5)  can  be 
written 

E(aq,B)  -  If(8-a,M)E(aqta)  *  E„(aq,8)  (C28) 

where  E(oq,a)  is  the  desired  piece  of  the  total  spectrum,  E „(ffq,0)  is  the 
rest  of  the  spectrum  the  influence  of  which  is  to  be  minimized,  and  H(0-a,A0) 
is  an  isolating  function.  The  isolating  function  has  a  value  of  1  in  a  region 
L8  wide  about  direction  a  and  is  zero  elsewhere.  Its  effect  is  to  allow  an 
artificial  distribution  of  the  energy  in  the  spike  of  E(cq,a)  over  a  finite 
arc  width  A 8  . 

19.  If  Equation  C28  is  substituted  into  the  integral  of  Equation  C27, 
one  obtains 
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Some  effort  is  required  to  evaluate  the  first  integral  on  the  right  side  of 
the  last  equality.  The  function  H(0-a,A$)  has  the  simple  effect  of  reset¬ 
ting  the  limits  on  the  integral  so  that  one  can  write 


f2"  £(*-«, AO  e *  •i*<7.-y„>l(rt 
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The  integrand  can  be  expanded  as  infinite  sums  of  the  products  of  Bessel 
functions  of  the  gage  separations  and  sinusoidal  functions  of  8  .  The 
sinusoids  can  be  integrated  and,  if  A 8  is  small  enough,  results  can  be 
expressed  as  products  of  A 8  and  sinusoids  of  a  .  The  infinite  sums  of  the 
resulting  products  of  Bessel  functions  of  the  gage  separations  and  sinusoidal 
functions  of  a  can  then  be  contracted  to  an  exponential  function. 

Equation  C30  then  becomes  the  approximate  expression 

e-ika[(JC,-xB)co«l  *  ^  Ik,  [<*.-*„>'<>•«  ♦  (7.-7„)eo».J  ^  (C31) 
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The  expression  for  the  frequency  cross -spectra  then  takes  the  form 
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20.  If  the  frequency  cross-spectra  of  Equation  C32  are  substituted  into 
the  model  Equation  C2S,  one  obtains 
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In  the  second  equality  of  Equation  C33,  E(oq,o )&$  was  moved  outside  the 
double  sum  because  it  does  not  explicitly  depend  on  m  or  n  .  In  the  third 
equality,  the  double  sum  and  weights  are  moved  inside  the  integral  because  the 
computation  can  be  done  in  any  order.  Also  in  the  third  equality,  the  weights 
and  the  exponential  functions  are  bracketed  together  with  the  double  summation 
signs,  thereby  grouping  the  summation  dependent  terms. 

21.  The  bracketed  terms  in  the  last  equality  of  Equation  C33  are  called 
window  functions ,  the  general  form  of  which  is  denoted  U(o ,8)  and  defined  by 

W(a,8)  -  £  £  <«*> 

n-X  b-1 

and  the  special  form  where  a  -  8  is 

»(«..)  -  £  £  ».(«)  -.‘(a)  <C35> 
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Using  the  window  function  notation,  the  estiaated  spectrum  from  the  last 
equality  of  Equation  C33  becomes 

E(at,a)AB  »  £(<rq.a)  A*  I/(a  ,a)  ♦  £*  W(a,  8  )  ^(aq,  9  )  d8  (C36) 

22.  A  constraint  on  the  problem  arises  where  there  is  no  "noise" 
spectrum  and  the  true  spectrum  is  all  that  exists.  That  is,  if 
£»(<*<!*«)  "  0  .  it  is  required  that  fi(oq,a) A0  -  E(<rq,a)&9  .  For  this  expres¬ 
sion  to  be  true,  W(a,a)  -  1  .  Note  that  the  expression  for  V(a,a)  can  be 
factored  as  follows 
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Because  the  term  in  the  second  bracket  of  the  last  equality  is  simply  the 
complex  conjugate  of  the  term  in  the  first  bracket,  Equation  C37  can  be 
written 
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and  it  is  seen  that  the  window  is  always  a  positive  and  real.  The  constraint 
that  U(a,a)  -  1  is  satisfied  (not  necessarily  uniquely)  if 

n  n 

£  wB(a)  eU<,<x"co“  *  7n‘i“)  -  £  wn*(a)  e-ik«(x"co“  *  .  1  (C39) 
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The  more  general  window  V(a,8) 

I. 

W(a,8)  - 


can  also  be  disassembled  and  rewritten  as 
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and  it  is  noted  that,  because  it  is  a  magnitude  squared,  its  values  are  real 
and  positive  or  zero  for  all  directions  8  . 

23.  Vith  V(a,a)  -  1  ,  our  estimated  spectrum  from  Equation  C36  becomes 

i(<7,.a)A*  =E(oq,a)&8  *  W(a,8)  E„(aq,8)  d6  D 

To  minimize  the  contribution  from  the  "noise”  spectrum,  the  integral  in 
Equation  C41  must  be  made  as  small  as  possible.  It  is  noted  that  E(<7q,a)  is 
of  fixed  value  and  positive,  that  V(a,8 )  is  positive  by  Equation  C40  and 
that  ER(crq,0)  is  positive,  so  that  everything  on  the  right  side  of  Equation 
C41  is  positive.  Making  the  integral  as  small  as  possible  is  then  equivalent 
to  making  £(crq,a)  as  small  as  possible.  That  can  be  done  for  any  expression 
defining  £(<7q,o)  .  Instead  of  minimizing  Equation  C41,  one  can  minimize 
Equation  C25,  which  is  a  more  basic  model  definition.  At  the  same  time  one 
must  satisfy  the  window  constraint  of  Equation  C39. 

24.  Minimization  means  finding  the  derivatives  of  Equations  C25  and  C39 
with  respect  to  the  real  and  imaginary  parts  of  the  complex  weights  wn  , 
setting  the  derivatives  equal  to  zero  and  solving  the  subsequent  set  of 
simultaneous  equations  for  the  weights.  Use  of  these  minimizing  weights  in 
the  model  Equation  C25  yields  a  maximum  likelihood  estimate  of  the  spectral 
density  at  frequency  oq  and  direction  a  .  Direction  binwidth  &8  contin¬ 
ues  to  be  undefined  until  the  minimization  has  been  completed. 

25.  Minimization  of  Equations  C25  and  C39  at  the  same  time  entails  use 
of  a  method  known  as  Lagrange  multipliers,  a  subject  covered  in  most  advanced 
calculus  texts  of  which  Hildebrand  (1965)  is  one.  To  set  it  up,  it  is  noted 
that  if  Equation  C39  is  true  then 

1  -  £  „n(a )  e*V*"eo“  *  .  o  (C42) 
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and  Equation  C42  can  be  added,  with  a  coefficient,  to  the  right  side  of 
Equation  C25  without  changing  the  truth  of  Equation  C25  or  Equation  C39.  This 
procedure  gives  us 
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The  coefficient  A  is  called  a  Lagrange  multiplier.  Upon  zeroing  the 
derivatives  of  Equation  C43  with  respect  to  the  real  and  imaginary  parts  of 
the  vB  ,  one  will  obtain  2N  equations  (  N  equations  for  the  real  parts  and 
N  for  the  imaginary  parts)  that  will  yield  the  wn  in  terms  of  A  .  If 
these  weights  are  used  in  Equation  C39,  another  equation  is  obtained  that  can 
be  used  to  solve  for  A  . 

26.  To  illustrate  this  procedure,  consider  an  array  of  3  non-collinear 
wave  gages.  The  analysis  will  work  with  any  larger  number  of  gages,  but  the 
bookkeeping  takes  longer  to  write.  By  showing  the  derivation  for  a  3-gage 
array,  one  can  see  the  generalization  for  a  larger  array.  With  3  gages,  there 
are  (  N  -  )  3  time  series ,  so  that  for  each  frequency  aq  ,  there  are  3 
Fourier  coefficients  and  three  weights  wn  to  find.  Each  weight  is  complex. 
The  real  and  imaginary  parts  of  the  weights  are  identified  by  defining 


and  its  complex  conjugate 
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Products  of  weights  as  they  appear  in  Equation  C43  take  the  form 

wn<  *  <wbR  +  iwol><w,*  ~ 
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To  save  space,  the  parentheses  and  the  independent  variables  they  carried  have 
been  dropped  for  now.  They  will  be  replaced  at  the  end  of  the  derivation. 

27.  The  six  differentiations  required  to  minimize  Equation  C43  can  be 
written  as 

,o)A*  -  -£-£(<7.0)4*  -  0  n  -  1,2,3  (C47) 

3*aa  ’  q 

or,  using  the  right  side  of  Equation  C43  and  expanding  the  differentials, 
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where  the  terms  in  the  curly  brackets  expand  to 
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+  (V2R  +  W2I  ^  E22  *  [  (w2BV38  *  W2IW3I^  +  ^  '•W2II,r3R  "*  W2RW3I^  1 E32 

*  [(W3RW18  *  W3IWll)  +  ^I^XR  '  W3RV1I>)E13 

*  [  (wznwzg.  *  W3IW2l)  *  *(W3IV2R  -  W3RW2l)]^23  +  (W3R  *  V31^E33 

♦  *[  1  -  ("a  ♦  iwu) e^11'""  *  r,«xon)  _  ^  +  *  V4-1 

-  (wJR  ♦  iw3I)eU«Cl>co-  * 


(C48) 


(C49) 
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Taking  the  partial  differential  of  Expression  C48  with  respect  to  wa  and 
setting  it  equal  to  zero  yields 


^VlsFll  *  ^V2SL~^V22^2X  *  (w3£-i*^3i)^3i  +  (Wa+^W2 1^X2 


(C50) 


Similarly  the  partials  with  respect  to  wa  and  vn  lead,  respectively,  to 


<WUt+iwlI>£21  +  <Wia~iwu)£U  *  2w2K£,22  *  <WaTiW3I>£32 


(C51) 


♦  (W2R+^W2I^£32  *  (wia~iwll)^13  (w2K~^V2I^£23 


*  2**  -  AeUt»l^eo“  *  T>ti“)  -  0 


(C52) 


The  partials  with  respect  to  wlx  ,  w2I  and  w3I  yield,  respectively, 


2wll£ll  *  ^W21*^W2s)^21  *  ^V31*^V3b)^3X  *  (v21~^v2i)^X2 


(C53) 


(vlZ~^wm)^21  +  (WU*iWlx)El2  *  2w2I£22  *  (V3I+^W 22^32 


(C54) 


♦  (Wjj-IWjjJEjj  -  iAetk,<x,C0,“  *  yj*ioa)  .  o 


(Wlj-iWj^JEjj  +  (w2I-iw2R)£32  ♦  (^♦IWjjjJEjj  ♦  (Vjji iW2*>E23 


(C55) 


2v3IE„  -  .  o 


Equations  CSO,  CS1  and  CS2  can  be  written  in  more  abbreviated  form  by  using 
the  definitions  of  Equations  C44  and  C45 .  If  terms  containing  X  are  also 
moved  to  the  right  sides  of  the  equations,  one  obtains 


S»«*u  *  »&>  *  •  -A,  *  *Ai  ■ 


(056) 


»Ai  *  ♦  “A.  ♦  »A,  ■  • ,-1”'  <C57> 
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(C58) 


VlEn  ♦  v:^2  +  VlEl 3  *  v2£23  ♦  2v3K£33  -  Ae*^'0*"  *  ,*,1“) 

If  Equation  CS3  is  multiplied  by  1  and  terns  are  rearranged,  the  following 
sequence  is  Induced 

2-*wi:£u  *  ^vtm)E2i  *  ^  (V3I+ ^V3s) EZ1 

*  i(w2I-lwffl)£u  «■  i(v3I-iw33l)£13  -  ia  A  ea»<x'eo“  *  y'#i“} 

“  ♦  (iw2I_w2S^£21  *  ^W3I~W3S^£31 

*  (i’'B~n>Sa  *  «»„«»)*„  *  - - 

■  2-*wn£u  -  (v2SL-lv2Z)E2l  -  (w3a-lw3I)£31 

*  (»a.i*„)£tt  *  <»,.!»„>£„  •  - - 


(C59) 


Writing  Equation  C59  in  abbreviated  form  (using  Equations  C44  and  C45)  and 
putting  the  term  containing  A  on  the  right  side  results  in 


2iv  E  -  v*E  -w'E  +vE  -  -A«U«<x,eo**  *  r',ina, 

zawllfill  *'2*'21  3C31  TU  *3^13  Ae 


(C60) 


The  same  set  of  operations  applied  to  Equations  C54  and  CS5  provides,  respec¬ 
tively, 


*A,  -  "i£u  *  -  •*»  -  "A,  -  •  ,1“”’ 

«A>  -  -  “,X,  - 


C  C61  > 
(C62) 


If  Equation  C60  is  subtracted  from  Equation  C56,  there  results 

<2*a  -  2itru)ru  -  2w*E21  *  2w3£3l  -  2Ae^Cx’eo“  *  <c63> 

or,  upon  putting  the  first  term  on  the  left  in  simplified  form  and  dividing 
by  2, 


Similarly,  subtracting  Equation  C61  from  Equation  CS7  leads  to 


(C64) 
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*;cu  ♦ 

and  Equation  C62  taken  from  Equation  C58  yields 

v*fu  ♦  v'Ejj  *  v3-f33  - 


(C65) 


(C66) 


Equations  C64,  C65  and  C66  are  very  similar  in  form  and  the  left  side  of  each 
can  be  written  as  a  sum.  On  doing  this,  the  equation  set  is  expressed 

3 


l >X,  - 


for  n  -  1,2,3 


(C67) 


If  one  multiplies  Equation  C64  by  wx  ,  Equation  C65  by  w2  and  Equation  C66 
by  w3  ,  there  is  obtained  a  new  set  of  equations  that  can  be  expressed  in  the 
summation  form  of  Equation  C67  as  the  product  of  both  sides  with  wn  ,  i.e., 


v  Y  w'E  •  w  Ae^<x"eo“  *  y»*i“) 
n  -  ^  a  an  n 


for  n  ■  1,2,3 


(C68) 


or,  upon  moving  wn  inside  the  sum  because  the  sum  is  only  over  index  m. 


for  n  ■  1,2,3 


(C69 ) 


If  all  three  equations  are  now  added  together,  a  very  simple  and  useful  result 
is  obtained.  Adding  the  equations  in  the  summation  expression  is  the  same  as 
taking  the  sum  over  index  n  of  both  sides  of  Equation  C69,  or 

3  3  3 


n-l  o-l 


Z  Z  vX  ■  Z  ^LVV°“ ' 

D«1 

3 

■  A  ]C  w  e1*9^60*"  *  r»*in“) 


(C70) 


The  double  sum  on  the  left  side  of  Equation  C70  is  recognized  as  the  estimate 
of  the  frequency-direction  spectrum  defined  in  Equation  C25.  The  sum  on  the 
right  side  of  Equation  C70  is  identified  as  satisfying  the  window  constraint 
expressed  by  Equation  C39  so  that  the  sum  can  be  replaced  by  1.  Equation  C70 
reduces  to  the  simple  expression 

E(,aq,a)  LB  •  \  <C71> 
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Hence,  if  A  con  be  found,  Che  estimator  will  be  complete.  Note  that  while 
the  number  of  gages  has  been  limited  to  3  to  keep  the  bookkeeping  simple,  all 
operations  so  far  have  been  linear.  Extension  of  this  derivation  to  more 
gages  is  straightforward,  consisting  mostly  of  changing  the  upper  limit  on  the 
summations.  This  derivation  will  continue  to  assume  3  gages,  again  for 
bookkeeping  simplicity,  but  the  result  is  more  general. 

28.  To  find  X  ,  it  is  noted  that  Equations  C64,  C65  and  C66  form  a  set 
of  simultaneous  equations  in  the  conjugates  of  the  weights  wn  .  These 
equations  can  be  written  in  matrix  form  as 


Eu  E21  E31 

vi 

eikq(x,co»a  ♦  y,*ina> 

E12  E22  E32 

• 

W2 

*  A 

^ikgCx^eoM  ♦  y3iina) 

&13  E23  E33 

J 

W3* 

glkq(x3coaa  ♦  y3*ina) 

v.  J 


(C72) 


To  find  the  weights,  both  sides  of  Equation  C72  must  be  multiplied  by  the 
inverse  of  the  matrix  Ej,,  .  This  inverse  is  computed  by  following  any  recipe 
for  matrix  inversion.  On  finding  the  inverse  and  doing  the  multiplication, 
Equation  C72  becomes 


•1 

W1 

C21 

E31 

• 

W2 

»  A 

Eu 

El 

Z22 

*3*2 

W3 

> 

*1*3 

W 

£,‘1 

C23 

> 

■  ' 

elk,(x,co»a  *  y,*lna) 

eika(xjeo»a  ♦  y3*ina) 
ii,(x3coia  *  yjiim) 


(C73) 


which  can  be  written  as  individual  equations  by  doing  the  matrix  multiplica¬ 
tion  to  form 


•  * 
vx  »  X 


^■IglkqtX^OM  ♦  y,»ina)  +  ^-lglkqCXjCOM  *  73«iix»>  +  ^-l^^CXjeo**  *  y^lna) 

XX  21  3X 


(C74) 


w2  *  X 


_-l  ik,(x,eo*a.y,«iaa)  _-l  tkqCx,eo»«*yj*iaa)  _-l  Utq(x3co««  .  y,iioa)  (C75) 
*12e  *  *22®  +  C.^26 


w3  -  A 


P-X  ikq(x,co*a  ♦  y,*inn)  _-l  Utq(x3eo»«  *  y3*ina>  _-X  Utq(Xjeo*«  ♦  y,iino) 
*13e  +  *23®  *  *33® 


(C76) 
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Equations  C74,  C75  and  C76  are  of  forms  involving  sums  of  similar  terms.  They 
can  be  compressed  into  a  single  summation  equation  that  looks  like 

vn*  -  for„  -  1,2,3  (C77) 

mml 

This  expression  is  readily  manipulated  to  isolate  X  .  If  both  sides  of 
Equation  C77  are  multiplied  by  *  7n*in“>  t  there  results 


3 

„  g-*V*«eo“  .  y„*ina>  ^  £  £-l  .  y^ina) 

•>1 


3 

^  yv  ^-leik,C*,so»«  *  y^iaa>e-li,(xnco«a  ♦  yn»ins> 
no 

«•! 


(C78) 


3 

A  p-lflika[(*B-*n)eoia  *  (y.-y„)«ina] 
jL*d  BD 


for  n  -  1,2,3 


where  the  exponential  multiplier  can  go  inside  the  summation  because  the  sum 
is  only  over  index  m  ,  and  the  exponents  are  combined  in  the  last  equality. 
If  the  3  equations  represented  by  Equation  C78  are  added  together,  i.e., 
summed  on  both  sides  over  index  n  ,  the  result  is 


3 

1 

n»l 


3  3 


£  j/*e'**Q^XnCO*°1  *  m  £'1eUt,,t(x"'*n,C0,a  *  lT»'T ",*inal 

n»l  »■ 1 

3  3 

W  £  ^-leli,rCx.-xn)co«a  *  (y.-yn)*ina] 


n-1  b>1 


(C79) 


The  left  side  of  Equation  C79  satisfies  the  window  constraint  of  Equation  C39 
so  that  it  equals  1  .  This  result  leaves 

3  3 

l  m  £  £-l  glJt,t(*.-xn)co»a  .  (y.-yn)iino)  (C80) 

n»l  sal 


or 
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3  3 

£  £  £-1  ei*g [<*,-*„)<»•«  *  <7.-7*>»ii»l 

■»!  »*l 


(C81) 


Hence,  an  expression  for  A  has  been  found,  and,  by  virtue  of  Equation  C71, 
so  has  an  estimate  of  the  frequency -direction  spectrum,  i.e., 


E(< 7q,a)  M 


_ 1 

E  E  «ik,t<1"‘x",co“ 

n-l  b-1 


(y«*ynJ,lnal 


(C82) 


This  is  the  result  being  sought,  a  maximum  likelihood  estimate  of  the  wave 
variance  per  frequency  increment  at  frequency  o,  and  look  direction  a  . 

29.  The  angle  increment  L8  still  has  not  been  defined,  however.  This 
definition  can  be  found,  symbolically  at  least,  by  noting  that,  in  practice, 
Equation  C82  is  commonly  solved  for  a  number  of  evenly  spaced,  discrete  lock 
angles  defined  by 

8t  •  r&6  for  r  »  1,2,_.,J?  (C83) 

such  that,  for  the  360  degrees  or  2jt  radians  of  arc  that  can  be  examined  with 
a  two-dimensional  array, 

■  RA8  »  2jt  (C84) 


from  which  M  -  2ir/R  .  With  this  definition  of  LS  and  the  set  of  discrete 
directions  given  by  Equation  C83,  the  frequency- direction  spectral  estimate  at 
frequency  and  direction  9t  is 


_ 1 _ 

2x  £  £  £-\  dlka((x.-xn)co**r  ♦ 

E  n»l  »-l 


(C85) 


Again,  nothing  has  been  done  in  this  derivation  that  precludes  adding  more 
gages  to  the  array.  Discussion  was  limited  to  three  gages,  mostly  so  the 
equation  sets  and  matrix  inversion  equations  would  fit  into  a  reasonable 
amount  of  space.  Extension  of  the  result  given  by  Equation  C85  is  achieved 
symbolically  by  simply  changing  the  upper  limits  on  the  sums,  so  a  more 
general  result  is 
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2»  £  £  ^,-1  .  (y.-y„).la»,J 

^  n»l  »il 


(C86) 


where  N  is  the  number  of  gages,  as  before. 

30.  One  muse  still  do  some  empirical  work  to  pin  down  AS  (or  R  ), 
because  the  directional  resolution  of  the  estimator  depends  very  strongly  on 
both  the  Array  design  and  the  nature  of  the  true  spectrum  that  governs  the 
data  collected  with  the  set  of  gages.  One  should  use  synthetic  data,  for 
which  the  true  spectrum  is  known,  to  examine  the  response  of  a  particular 
array  design.  From  such  results,  one  can  determine  the  extent  to  which  the 
array  can  resolve  various  wave  approach  directions,  including  cases  where 
waves  come  from  several  directions  at  the  same  frequency.  Some  very  good 
guidance  on  proper  placement  of  gages  in  array  design  is  given  by  Davis  and 
Regier  (1977). 

31.  Once  an  array  is  designed  and  deployed  such  that  data  are  obtained, 
one  applies  this  method  by  computing  the  cross-spectra  between  all  gage  pairs, 
forming  the  cross -spectral  matrix  E,*,  for  each  frequency  <7q  ,  computing  the 
inverse  of  the  cross-spectral  matrix,  and  applying  the  estimator  Equation  C86 
for  each  discrete  direction  0r  . 

32.  This  method  is  not  constrained  to  conserve  total  variance,  so  when 
the  frequency-direction  spectral  estimates  are  integrated  with  respect  to 
direction,  the  frequency  spectrum  is  not  necessarily  recovered  completely. 
There  are  two  common  ways  to  deal  with  this.  The  simplest  one  is  to  perform 
the  integration  with  respect  to  direction  at  frequency  o,  ,  compute  the  ratio 
of  the  variance  from  the  frequency  spectrum  to  the  result  of  the  integration, 
and  multiply  all  values  of  the  frequency-direction  spectrum  at  this  frequency 
by  the  computed  ratio.  This  form  of  correction  distributes  the  error  as  a 
uniform  percentage  over  the  entire  directional  distribution.  A  second 
correction  method,  known  as  Iterativ  .  Maximum  Likelihood  Estimation,  is  based 
on  a  method  derived  by  Pawka  (1983).  The  reasoning  of  Pawka's  method  is  that 
if  the  MLE  result  is  fed  back  into  the  estimation  algorithm,  the  result  of  the 
second  pass  tends  to  be  an  attenuated  form  of  the  first  pass.  If  results  from 
the  first  pass  are  multiplied,  direction-by-direction,  by  the  ratio  of  first 
pass  results  to  second  pass  results,  the  new  estimate  tends  to  be  closer  to 
the  true  spectrum  than  the  first  pass  results.  Extensive  empirical  tests  by 
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Pawka  indicate  that  repeated  application  of  this  procedure  tends  to  produce  an 
estimator  that  is  much  better  than  the  single-pass  MLE,  not  only  giving  better 
fidelity  to  the  total  variance,  but  also  sharpening  resolution  in  the  detailed 
structure  of  the  directional  distribution  of  variance.  The  mathematics  of 
this  procedure  are  not  discussed  here  because  they  are  rather  well  described 
by  Pawka  (1983)  and  by  Oltman-Shay  and  Guza  (1984). 

33.  Thus,  this  appendix  has  outlined  some  of  the  reasoning  and  some  of 
the  mathematical  detail  of  maximum  likelihood  estimation  of  wave 
directionality  when  used  with  a  two-dimensional  array.  Equation  C86  was  used 
identically  to  estimate  the  frequency-direction  spectra  discussed  in  the  body 
of  this  report. 
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Appendix  D :  Notation 


[  ]*  Complex  conjugate  of  [  ] 

a  Wave  amplitude 

a,,  Amplitude  of  11th  component  wave 

A  Duration  of  time 

Aq  Fourier  series  cosine  coefficient  at  q1*  frequency 

B  Span  along  x- direction 

Bq  Fourier  series  sine  coefficient  at  q*11  frequency 

C  Span  along  y-direction 

C,,,,  Coincident  spectral  density  between  gages  m  and  n 

d  Water  depth 

dd  Mnemonic  for  two-digit  day  of  the  month 

D  Complex  wave  amplitude 

D(fa,0a)  Element  of  directional  distribution  function  at  XXth 

frequency  and  111th  direction 

Da  Complex  amplitude  of  n**  component  wave 

E  Frequency- direction  spectral  density 

£  Estimate  of  frequency-direction  spectral  density 

E^  Cross-spectral  density  between  gages  m  and  n 

C  Element  of  row  m  and  column  n  of  matrix  inverse  of 

E- 

E,  Noise  component  of  frequency- direction  spectrum 

f  Cyclic  frequency 

fn  Cyclic  frequency  of  n*11  component  wave 

fp  Spectral  peak  frequency 

f p Peak  frequency  of  incident  spectrum 
f r  Peak  frequency  of  reflected  spectrum 
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s 

H 

H-.i 

H»,r 

k 

m 

mm 


M 

n 

N 

q 

Q 

Q» 

r 

R 

S 

S(f) 
S(fn) 
S(f a,8m) 

S(«) 


ss(fa) 

t 


Gravitational  acceleration 
Isolating  function 

Spectrum-based  characteristic  wave  height 
Characteristic  height  of  incident  spectrum 
Characteristic  height  of  reflected  spectrum 
Radian  wavenumber 
Wavenumber  of  n^  component  wave 
Index  integer 

Mnemonic  for  two-digit  representation  of  month  in  a 
calender  year 

Upper  limit  of  indices  denoted  m 
Index  integer 

Upper  limit  of  indices  denoted  n 
Index  integer 

Quadrature  spectral  density 

Quadrature  spectral  density  between  gages  m  and  n 
Index  integer 

Upper  limit  of  Indices  denoted  r 
Covariance  function 

Spectral  density 

Frequency  spectral  density 

Frequency  spectral  density  at  n6*1  discrete  frequency 

Frequency -direct ion  spectral  density  at  n61*  discrete 
frequency  and  111th  discrete  direction 

Direction  spectral  density 

Incident  frequency  spectral  density  at  rf*  discrete 
frequency 

Reflected  frequency  spectral  density  at  n^  discrete 
frequency 

Time 


DA 


T  Arbitrary  displacement  In  time 

Tp  Spectral  peak  period 

Tp(1  Peak  period  of  incident  apectrua 

Tpit  Peak  period  of  reflected  apectrua 

vm  Complex  weighting  function  for  Fourier  coefficient 

contribution  from  gage  a 

w^r  Imaginary  part  of  va 

Real  part  of  va 

wa  Complex  weighting  function  for  Fourier  coefficient 

contribution  from  gage  n 

w^j  Imaginary  part  of  wn 

Wpjt  Real  part  of  wB 

U  Window  function 

x  Horizontal  cross* shore  coordinate 

x.  Cross* shore  coordinate  of  gage  a 

x,j  Cross -shore  coordinate  of  gage  n 

X  Displacement  or  lag  in  x- direction 

Xq  Element  of  discrete  Fourier  transform  at  q11*  frequency 

y  Horizontal  longshore  coordinate 

yy  Mnemonic  for  two-digit  representation  of  year 

yymadd  Date  mnemonic,  a  concatenation  of  year,  month,  and  day 

mnemonics 

ya  Longshore  coordinate  of  gage  m 

ya  Longshore  coordinate  of  gage  n 

Y  Displacement  or  lag  in  y-direction 

a  Wave  direction  as  fitted  by  maximum  likelihood 

estimation 

At  Increment  or  step  in  time 

bo  Increment  of  radian  frequency 

b$  Increment  of  wave  direction 


D5 


c  Initial  wave  phase 

«B  Initial  phase  of  11th  component  wave 

1 j  Sea  surface  displacement 

Ht  Element  of  sea  surface  displacement  time  series  at 

time  step 

X  Wavelength 

Lagrange  multiplier 

8  Wave  direction 

8a  Element  m  of  a  discrete  sec  of  wave  directions 

9a  Direction  of  n**  component  wave 

8 _  Peak  wave  direction  - - — %  rw_ 

0pi  Peak  direction  of  incident  wave  spectrum  J _ 

8p  t  Peak  direction  of  reflected  wave  spectrum 

a  Radian  wave  frequency 

oa  Radian  frequency  of  n**  component  wave 

oq  Element  q  of  a  set  of  discrete  radian  wave  frequencies 

X  Reflection  coefficient  based  on  wave  heights 

Xao  Reflection  coefficient  based  on  characteristic  wave 

heights 

Xp,i  Reflection  coefficient  at  the  peak  frequency  of  the 

incident  wave  spectrum 

Xp>r  Reflection  coefficient  at  the  peak  frequency  of  the 

reflected  wave  spectrum 
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